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The universal coupling of matter and gravity is one of the most important features of general 
relativity. In quantum gravity, in particular spin foams, matter couplings have been defined in the 
past, yet the mutual dynamics, in particular if matter and gravity are strongly coupled, are hardly 
explored, which is related to the definition of both matter and gravitational degrees of freedom on 
the discretisation. However extracting this mutual dynamics is crucial in testing the viability of the 
spin foam approach and also establishing connections to other discrete approaches such as lattice 
gauge theories. 

Therefore, we introduce a simple 2D toy model for Yang-Mills coupled to spin foams, namely 
an Ising model coupled to so-called intertwiner models defined for SU(2)fe. The two systems are 
coupled by choosing the Ising coupling constant to depend on spin labels of the background, as 
these are interpreted as the edge lengths of the discretisation. We coarse grain this toy model via 
tensor network renormalization and uncover an interesting dynamics: the Ising phase transition 
temperature turns out to be sensitive to the background configurations and conversely, the Ising 
model can induce phase transitions in the background. Moreover, we observe a strong coupling of 
both systems if close to both phase transitions. 

PACS numbers: 04.60.Pp,04.60.Nc,05.50.+q 

I. INTRODUCTION 

The universal coupling to all types of (known) matter 
is one of the most profound features of general relativity. 

This property has proven to be invaluable in testing the 
theory, make progress in cosmology and interpret phe¬ 
nomena such as black holes as originating from the grav¬ 
itational collapse of matter. However, general relativity 
as a classical and thus deterministic theory of a dynam¬ 
ical spacetime is at odds with the fundamental quan¬ 
tum nature of matter, which is described by quantum 
field theories on a fixed, not necessarily flat, background 
spacetime. 

It is expected that this inconsistency can be cured by 
quantising gravity [1]. Inevitably this not only requires 
this theory of quantum gravity to possess a classical limit 
consistent with vacuum general relativity, but it must 
also be possible to couple matter to it and obtain a dy¬ 
namics for the composite system. Most importantly, this 
should allow us to study systems, for which gravity and 
matter are strongly coupled, e.g. close to the Big Bang, 
where one can expect to find new physics. However, these 
regions are currently not accessible by experiments. On 
the other hand, the theory must not conflict already well- 
tested physics. Hence, it should also be possible to iden¬ 
tify a regime, in which the effective dynamics (of quan¬ 
tum matter and gravity) behave like quantum field theory 
on an almost fixed background spacetime, e.g. consistent 
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with data observed at accelerators. Moreover, the cou¬ 
pling of matter might reveal new properties and deepen 
our understanding of the theory of quantum gravity, as 
it has been demonstrated e.g. in [2] in the context of 
asymptotic safety [3]. 

The status of coupling matter to candidate theories 
of quantum gravity is as diverse as the theories them¬ 
selves. For approaches relying heavily on quantum field 
theory methods, e.g. asymptotic safety [3], in particu¬ 
lar choosing (albeit not necessarily hxing) a background 
spacetime, coupling matter to gravity is (at least con- 
ceptionally) straightforward, whereas theories relying on 
(either auxiliary or fundamental) discretisations, such as 
spin foam models [4, 5], causal dynamical triangulations 
[6] or causal sets [7], additionally face the issue of defining 
gravity and matter on the same discrete structure. 

A candidate theory of quantum gravity defined on a 
discretisation, so-called spin foam models [4, 5], moti¬ 
vate this paper. Spin foams are a path integral approach 
related to loop quantum gravity [8, 9] and are defined on 
a 2-complex, i.e. a collection of vertices, edges and faces, 
which is frequently chosen to be dual to a triangulation. 
These elements of the foam carry representation theo¬ 
retic objects of the underlying symmetry group; for Eu¬ 
clidean signature this is SO(4), whereas SL(2,C) is cho¬ 
sen for Lorentzian signature. Irreducible representations 
can be found on the faces, whereas on the edges projec¬ 
tors onto an invariant subspace (in the tensor product of 
vectorspaces of representations meeting at this edge) are 
located. The spin foam model assigns an amplitude to a 
configuration of the foam, i.e. a colouring with represen¬ 
tations and intertwiners, where one obtains the full am- 
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plitude after summing over all configurations. If the 2- 
complex has a boundary, the associated boundary Hilbert 
space is that of loop quantum gravity on the boundary 
graph^; thus a spin foam assigns an amplitude to one 
possible transition of loop quantum gravity states. 

Readers familiar with lattice gauge theories may recog¬ 
nize that the data carried by a spin foam is very similar to 
the strong coupling expansion of lattice gauge theories. 
Indeed, spin foam models can be understood as gener¬ 
alized lattice gauge theories, since they carry projectors 
onto a smaller invariant subspace instead of the full Haar 
projectors. In the usual holonomy representation, this re¬ 
quires the introduction of additional group elements as¬ 
signed to the faces of the discretisation [14, 15]. In fact, 
these similarities can already be noticed for the classical 
continuum action, where Yang-Mills and gravity can be 
written as topological BF theory [16] plus constraints, 
which break the too many symmetries of the topolog¬ 
ical theory. The crucial difference is however that the 
simplicity constraints^ in the Plebanski formulation of 
general relativity [17] are manifestly background inde¬ 
pendent, whereas for BF-Yang-Mills [18] the constraints 
explicitly depend on a background metric. 

The before mentioned similar mathematical structure 
makes pure lattice gauge theories a very natural choice 
for (tentatively) examining matter couplings to spin 
foams, where one can focus on the conceptual issues of 
the coupling of matter and gravity, which is another mo¬ 
tivation for the present paper. In fact several different 
ways to couple Yang-Mills to spin foams have already 
been defined in the literature: In [19-21] the system of 
gravity and matter is described by a topological quan¬ 
tum field theory, in which propagating degrees of free¬ 
dom arise either by symmetry breaking in a low energy 
phase or matter is represented by defects in the geome¬ 
try. Conversely, [22, 23] define lattice gauge theory on a 
triangulation (dual to the foam) by modifying the Wilson 
action to depend on the geometric data provided by the 
foam, which can be interpreted as a geometry dependent 
coupling constant - an idea we will revive and test in 
this work for a toy model. In contrast to this, [24] de¬ 
rived a discrete action for (BF-) Yang-Mills [18] coupled 
to topological 3D gravity from the classical action and 
pointed out that in order to consistently make contact 
with boundary Hilbert spaces of the canonical theory, 
both theories ought to be discretized on the same dis- 


^ As most spin foam models are defined for triangulations, this 
does not contain all boundary graphs allowed in loop quantum 
gravity. See [10] for a spin foam model resolving this issue and 
also [11] on the challenge of linking loop quantum gravity and 
spin foam models. On the other hand, in [12, 13] a new rep¬ 
resentation of loop quantum gravity based on graphs dual to 
triangulations has been defined, which could be closer to spin 
foam models. 

^ Implementing these constraints in spin foam models results in the 
projectors onto smaller invariant subspaces, which is however not 
uniquely defined. 


Crete structure, here the 2-complex. In yet a different ap¬ 
proach [25] , Yang-Mills and gravity are unified in a larger 
symmetry group, which possesses SO(4) (or SL(2,C)) as 
a subgroup and the remainder of the group is associated 
to the Yang-Mills degrees of freedom. In the context of 
coupling fermions to spin foams [26], the introduction of 
U(l) Yang-Mills has been discussed in [27]. However, de¬ 
spite this multitude of interesting ideas, it has not been 
possible to thoroughly study the dynamics of the coupled 
system or even establish a link to usual lattice gauge the¬ 
ories defined on a flat background lattice. 

Similar to the pure gravity case, this issue is rooted in 
the introduction of the discretisation: Both for gravity 
and matter, the theory is not uniquely defined, e.g. for 
Yang-Mills one can choose the Wilson action or the heat 
kernel action [28], but moreover the crucial diffeomor- 
phism symmetry of general relativity is broken [29-33] 
and the theory depends on the chosen discretisation. In 
fact this also includes the chosen coupling (mechanism) 
of the two theories. However, physical predictions such as 
expectation values of observables, should not depend on 
the chosen regulator, which inevitably raises the ques¬ 
tion how to remove it. In spin foams alone this is a 
topic of debate [34], where two main options have been 
advertised: the first considers a sum over all discreti¬ 
sations (and even topologies as put forward by group 
field theories [35]), which (at least formally) removes the 
cut-off, whereas the second one proposes a refinement- / 
coarse-graining procedure that relates theories, bound¬ 
ary states and observables defined on different discreti¬ 
sations and can be understood as a (generalized) renor¬ 
malization procedure [36-46], since a notion of ‘scale’ is 
absent^. In this work we follow the latter idea, which 
eventually requires us to renormalize the coupled sys¬ 
tem, e.g. the EPRL-model [48, 49] coupled to U(l)- 
or SU(3)-Yang-Mills theory, which clearly is beyond the 
scope of this paper and beyond the possibilities of cur¬ 
rent renormalization techniques, despite recent progress 
in coarse graining algorithms for lattice gauge theories 
[50] and analogue spin foam models [51-54]. 

Instead, the goal of this paper is to introduce a toy 
model, in fact similar in spirit to the work [55, 56] for 
2D causal dynamical triangulations, that shares several 
crucial properties with spin foams and lattice gauge the¬ 
ories, but is considerably simpler such that we can study 
its renormalization group flow, here via tensor network 
renormalization [57, 58]. As the physical system we are 
interested in, our toy model consists of two parts: a ‘spin- 
foam inspired’ background, given by ‘intertwiner mod¬ 
els’ [59], which are based on the quantum group SU(2)fc 
[60, 61], and a Z 2 Ising model chosen for the matter part. 


® In contrast to standard lattice gauge theories, where the lattice 
spacing a is inherited from the background spacetime, geometric 
information has to be inferred from the geometric variables of 
the spin foam. One already faces a similar issue if one considers 
irregular lattices [47]. 
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Both systems can be understood as a dimensional reduc¬ 
tion of their 4D counterparts, in which the local symme¬ 
try is replaced by a global symmetry. Thus both systems 
pertain a crucial ingredient, namely intertwiners, which 
are located on the vertices of the lattice"'. We choose 
this lattice to be a regular square lattice. Then, these 
two systems are coupled by choosing the Ising coupling 
constant /? to be a function of background spin labels 
j®, similar to [22, 23], based on the interpretation that 
the background spins give the distance between two Ising 
spins. This idea actually resonates with the idea of local 
coupling constants for lattice gauge theories defined on 
irregular lattices in [47] , where the coupling constants are 
assigned to parts of the discretisation, e.g. edges or faces. 
Therefore we study the dynamics of coupled intertwiner 
degrees of freedom and are particularly interested in how 
the two systems affect each other. To do so, we study 
the Ising model on different backgrounds (including su¬ 
perpositions) for two different ways of coupling the two 
systems and extract phase diagrams of the model. 

The paper is organized as follows: in section II we in¬ 
troduce both the background and the Ising model and 
two different ways to couple the two theories. Section 
III describes the coarse graining algorithm. We split the 
results in two parts: section IV concerns the Ising model 
(for different temperatures) on a topological background, 
i.e. one fixed point intertwiner for the background, in 
section V we superpose several fixed point intertwiners. 
We close with a summary and discussion in section VI. 
In appendix A we have included some additional infor¬ 
mation on the quantum group SU(2)fe and the graphical 
calculus developed in [54, 59] helpful to understand some 
calculations in this article. Appendix B explains tensor 
network renormalization in more detail and gives a more 
thorough derivation of the formulae. 


A. The ‘gravitational’ part — Intertwiner models 


On the ‘gravity’ side, we pick the so-called intertwiner 
models defined in [59] by Dittrich and Kaminski as our 
dynamical background. These models are defined on a 
graph with 3-valent vertices with an underlying symme¬ 
try governed by the quantum group SU(2)fc at root of 
unity, i.e. the quantum deformation of the universal en¬ 
veloping algebra Uq{su{2)) with deformation parameter 
q = exp see e.g. [60, 61] and appendix A for more 
details, fc e N is also called the level of the quantum 
group and dehnes a natural upper cut-off on the spins 
j with jmax = Each 3-valent vertex is dual to a 
triangle, such that each of its links pierces exactly one 
edge of a triangle. Each of the links, and hence also the 
dual edges, carries a spin j £ ^ and a magnetic index 
m £ ^ with —j < m < j . To the vertex itself, we 
assign an amplitude: The dependence on the magnetic 
indices and triangle inequalities for the spins {je} are en¬ 
coded in the Clebsch-Gordan coefficients of 

the quantum group SU(2)fc. Additionally, we assign a 
factor a{ji, j 2 , ja) to the vertex, which only depends on 
the representation labels. Thanks to the imposed trian¬ 
gle inequalities, we will interpret the spin labels as the 
length of the edges of the triangulation. 

To sum up, the basic building blocks of our model - 
expressed in a graphical notation® - are: 


^(jl V2 J Js) 7712 m 3 





ji h 



( 2 . 1 ) 


II. THE TOY MODEL - ISING SPINS ON A 
DYNAMICAL BACKGROUND 

In this section we introduce the toy model of interest: 
in section IIA we Hrst present the ‘gravitational’ part fol¬ 
lowed by the (uncoupled) Ising model in section IIB. In 
section IIC we discuss the coupling of the two systems 
in more detail, including assumptions and a brief discus¬ 
sion on choosing the dual discretisation related to [22] 
and [24]. More importantly, we introduce two different 
couplings in sections IICl and IIC2, for which we can 
already deduce a very different qualitative behaviour. 


32 jl 32 jl 

= a{ji,j2,ja) 

33 33 


( 2 . 2 ) 


where q denotes the complex conjugate of the deforma¬ 
tion parameter q. 

In general, the factors a{{je}) and a'{{je}) can be cho¬ 
sen freely, but we are rather interested in the particular 
case of triangulation independent, i.e. topological, mod¬ 
els. For a model to be triangulation independent, the pre¬ 
dictions of the theory, e.g. expectation values of observ¬ 
ables and the partition function, must be independent on 


As an Abelian theory the Ising model only has trivial intertwiner 
spaces, an extension to non-Abelian (finite) groups is possible, 
see e.g. [62]. Note also that the intertwiner models [59] are also 
simpler than spin net models [51, 53, 54], which are also called 
analogue spin foam models. 

^ This possibility has been pointed out to us by Bianca Dittrich, 
who also raised the idea to examine a coupling ^{je) 


The peculiarities of the quantum group require us to specify a 
preferred direction, here the edges of the graphs come with an 
orientation pointing from bottom up. This is not crucial for 
most part of this work, thus we have included more details in 
appendix A. See also [54] and [59] for more thorough discussions 
of this graphical notation. 
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the chosen triangulation of the manifold. As shown by 
Pachner [63, 64], two triangulations of the same manifold 
are related to one another by a consecutive application 
of certain local changes of the triangulation, known as 
Pachner moves. 

There exist three Pachner moves in 2D: the 2-2 Pach¬ 
ner move, the 3-1 Pachner move and its inverse, the 1-3 
Pachner move. These moves are illustrated e.g. in equa¬ 
tions (2.5) and (2.8). The 2-2 move essentially exchanges 
the edge shared between two triangles by an edge con¬ 
necting the ‘tips’ of the original triangles. This gives two 
new triangles, such that the total number of 2-simplices 
in the triangulation is unchanged. For the 3-1 move on 


the other hand, one considers three triangles, where all 
three triangles share a common vertex. These three tri¬ 
angles are combined into one triangle by removing this 
(interior) vertex resulting in a coarser triangle; the in¬ 
verse 1-3 move reverts the process by placing a vertex in 
the center of the triangle and connecting it to the vertices 
of the coarse triangle. 

In the systems investigated in this work, invariance un¬ 
der Pachner moves directly translates into conditions on 
the amplitudes associated to the building blocks / trian¬ 
gles, in particular a({je}). Using the notation introduced 
above, we can write these conditions pictorially (see [59] 
for more details): 


• ‘Tilting’ condition: 



js 


which leads to the following relation for 



•^4 J3 ji 


the factors a({j}): 

^623 ®641 = 0362 0^54 


• 2—2 Pachner move: 


E 




Together with the following identity on the Clebsch-Gordan coefficients [59] 




(2.3) 


(2.4) 


(2.5) 


( 2 . 6 ) 


where the symbol in the square brackets is a modified {6j} symbol and dj = [2j + 1] denotes the quantum 
dimension of spin j, which is the quantum number of the classical dimension. Again, see appendix A for more 
details. One obtains the following condition on the factors a({j}). 



42 

44 



• 3—1 Pachner move: 



(2.7) 


( 2 . 8 ) 
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which gives the following equation for the factors a({j}): 


^23 — ^ 


Jl J2 J3 

k k J6 


( — 06510624^543 


and 


E 

i4,j5 J6 




which gives the following equation for the factors a{{j}): 


0 0123 = E 

44 , 45,46 


jl j2 js 
ji k k 


(_ 1 ) 44 + 46-43 

x/dja dj^ 


0641 0625 0453 


(2.9) 


( 2 . 10 ) 


( 2 . 11 ) 


Even if the factors a{{je}) (and a'{{je})) satisfy these 
conditions, i.e. the theory is topological, these factors 
are not uniquely determined. Further input is required, 
which usually is a restriction on the allowed spins. In [59] 
a wide variety of these models has been identified, here 
we briefly review the ones most interesting in this work. 
Note that from now on, we restrict ourselves to integer 
spins, j S N, i.e. representations of 80(3)^. Moreover, 
for all these topological theories, a({je}) = a'{{je}), such 
that we can simplify the notation: 

• Maximal spin J < jmax- 

This is the class of topological models we will be 
mostly interested in in this work. They are labelled 
by a spin J S N, which indicates that all spins 
j < J are allowed and excited, that is all couplings, 
which are allowed by the quantum group, involving 
these spins are allowed. J can be smaller or equal to 
the maximal spin of the quantum group, jmax = f. 
The corresponding factor a'^({je}) is given by 


dimension = 1. Hence jmax can only cou¬ 

ple with itself to the trivial representation j = 0; 
jmax ® jmax = 0. One Can construct a topological 
theory where only these two representations are al¬ 
lowed by choosing a({j}) as follows [59]: 

o(0,0,0) =a'(0,0,0) = 1 

^(jmaxj jmaxj 0) (jmax, jmax; 0) ~ 1 . (2.13) 

• Only j = 0, j = and j = jmax excited: 

For this model, we consider A: to be a multiple of 
four, k = Al. Then jmax = 2/ and = 1. As 
in the previous example, jmax = 2^ is of quantum 
dimension one, thus again all factors a({j}) which 
just contain j = 0 and jmax are equal to one (unless 
the coupling is forbidden), see (2.13). Thus the only 
remaining factors are those containing j = l\ Those 
factors a({j}) containing either j = 0 or j = 21 are 
proportional to the quantum dimension ol j = 1: 


a\ji,j2,k) = yEiE^v^EE^v^EiE^ x 


1 [2ji + l][2j2 + 1] 

jl j2 k 

/ [■/+!] 

J J J 

.2 2 2 - 


( 2 . 12 ) 


Hence, by choosing J, we have a direct handle on 
the maximal spin allowed on the edges of the lat¬ 
tice. Moreover, it is straightforward to consider 
linear combinations of these theories as it has also 
been studied in [54], albeit for more complex spin 
net models. 

• Only j = 0 and j = jmax excited: 

For even levels k of the quantum group, the maxi¬ 
mal spin jmax = I is an integer and is of quantum 


a(0, 1,1) = a{l, 0, 1) = a{l, 1, 0) = \/di 

a{2l,l,l) = = Vdi ■ (2.14) 

Thus the only non-trivial amplitude remaining is 
a{l, 1,1), which is given by: 

= a'(1,1,1) = \/di{di - 2) . (2.15) 

Note that these types of fixed points cannot be con¬ 
structed for all k = Al, since in a few cases the 2-2 
move cannot be solved, e.g. for A: = 12, A: = 24 or 
k = 32. See [59] for more details. Therefore we re¬ 
strict ourselves to the working cases, namely A; = 4, 
which is identical (up to signs) to the model J = 2 
discussed previously and k = 8. 
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1. Intertwiner models for f-valent vertices: 

The discussion so far only concerned 3-valent inter¬ 
twiners. In order to define the background for a square 
lattice, we have to define 4-valent intertwiners. 

In general, one can construct 7V-valent intertwiners 
from 3-valent ones by choosing a 3-valent graph with 
N outer edges and specifying additional N — 3 internal 
representations^. In our case, the 3-valent vertices come 
with additional restrictions on the allowed irreps, which 
we have encoded in the factors a({j}). These restric¬ 
tions translate into restrictions on the 7V-valent inter¬ 
twiners, if only the allowed representations for 3-valent 
intertwiners appear in the graph. Thus the set of allowed 
representations for 3-valent vertices determine those for 
A'^-valent ones; this is known as Reisenberger’s construc¬ 
tion principle [65], see also [54]. Note that different re¬ 
coupling schemes generically result in different N^-valent 
intertwiners. 

In the 4-valent case, there exist three different recou¬ 
pling schemes: 



(2.16) 

Due to our choice of factors a({j}), we see that the re¬ 
coupling schemes A and B are related to one another by 
a 2-2 Pachner move and thus give the same set of al¬ 
lowed representations and intertwiners. The case C, as 
it is also discussed in [54], is not satished for all models, 
in particular not the ones labelled by J, due to braiding 
in the quantum group, but still gives the same set of al¬ 
lowed representations as cases A and B. Nevertheless, in 
this work we will restrict ourselves to the planar case, i.e. 
cases A and B. 

Crucially, any choice of factors a({j}) giving a trian¬ 
gulation independent model, specihes a fixed point of the 
renormalization scheme (if only the background is consid¬ 
ered) we will use in this work, and thus a special family 
of intertwiners, |i^) in the 4-valent case. Given this, we 
define the projector on the 4-valent vertices, as: 

n"=ka)(4| , (2.17) 

which satisfy the projector condition V ■ V = V hy defi¬ 
nition. We can also define superpositions of intertwiners, 
which will however generically not satisfy this projector 
condition any more: 

= , (2.18) 


^ The graph is also called a recoupling scheme. Additionally we 
assume multiplicity free representation categories, i.e. in a tensor 
product of two irreps, each irrep appears at most once. 


where {cti} are the coefficients of the linear superposi¬ 
tion. To put this in context to the models under discus¬ 
sion here, let us express V for any choice of triangulation 
independent factors a({j}): 

h ji 

1 V 

T^vi{j}Am})=J2^ j, 

" " A 

js k 



J 3 J4 


To check the projector property for (2.19), the projectors 
have to be concatenated vertically (due to the special 
direction of the quantum group), see also [40, 59] for more 
details. For the rest of the work, we drop the superscript, 
since we always work on a square lattice. Eventually, we 
have the following partition function for the quantum 
gravity inspired background: 

Z= T’ii({je}eDii, {WejeD?;) ■ (2.20) 

B. The ‘matter’ part — Ising spins 

The goal of this work is to couple a matter system to 
this gravitationally inspired, possibly dynamical, back¬ 
ground. The simplest, non-trivial option ‘living’ on a 
lattice is the Ising model, a Z 2 spin system. On a given 
lattice, here with equidistant vertices, it is defined as fol¬ 
lows: 

On the vertices v of the lattice, one places group el¬ 
ements G Z 2 ; their interactions, restricted to nearest 
neighbours, are described by edge weights assigned to the 
edges e: 

^e{9,(e)9i{l)A) = 9s(e)9i(l)} : (2-21) 

where the inverse ‘temperature’ /3 ~ ^ is the coupling 
constant of the theory, G {—1,1} and the group 
multiplication is simply given by scalar multiplication. 
Each edge e comes with an orientation, where s(e) and 
t{e) denote the source and target of the edge e respec¬ 
tively. Inverting the orientation of an edge results in in¬ 
verting the argument of the corresponding edge weight, 
<^ 6 - 1 ( 5 ) = ^e(ff~^). Since 5 = —1 is its own inverse, 
changes of orientation do not affect the Ising model. The 
introduction of an orientation is rather a choice of nota- 
tional convenience and necessary for the quantum groups. 

Given these ingredients, the Ising model is defined via 
its partition function 

^=En exp{/3 9s(e)9t(l)} : (2.22) 

ISv} e 
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where the sum is over all configurations {gv} of group 
elements assigned to the vertices. This sum of configu¬ 
rations is weighted with the edge weights, which assign 
a higher weight to aligned spins, i.e. ( 7 s(e) = gt(e)- (2.22) 
also has a global symmetry: if all spins are simultane¬ 
ously flipped, i.e. = ±1 —>■ = TlVw, the partition 
function is unchanged. 

As it is well-known from Onsager’s solution [66], this 
model exhibits two different phases, separated by a sec¬ 
ond order phase transition at the critical inverse tem¬ 
perature /3crit ~ 0.4406..., in fact two different effects 
are dominating in the two phases. For f3 > /3crit, be. 
low temperature, the system is in the ordered phase, in 
which neighbouring spins are more likely to be aligned to 
minimize the energy of the configuration. For /? < /3critj 
i.e. high temperature, the spins essentially decouple and 
are not aligned any more. Still energy minimization and 
hence spin alignment is preferable, yet the higher the 
temperature, the smaller the energy gap becomes. Since 
there are many more configurations for which the spins 
are not aligned, they dominate the partition function. 

For the purpose of this work, it will be useful to briefly 
discuss the (group) Fourier transform of the Ising model, 
see also [51, 62]. Any function u}{g) on Z 2 can be ex¬ 
panded in characters y®: 


{9) = ^ ^k Xk{g) , 

k 

(2.23) 

^{g) Xkig) ■ 

(2.24) 


3 


k S {0,1} denote the irreducible representations of Z 2 
and Xk{g) is the character of the group element g in the 
representation k. Xk denotes the complex conjugate char¬ 
acter. This character is given by: 


Fourier components: 


Iinexp{/3'?s(e)<7t(e^} 

Isv} e 

( 2 . 28 ) 

^ Xfce( 5 s(e) 5 t(e)) 

{sAffcel e 

( 2 . 29 ) 

{fc.} e {£,„} V eDv 

( 2 . 30 ) 

{fee} e \eDv / 

, ( 2 . 31 ) 


where E denotes the number of edges in the lattice. We 
have used the fact that Xkig • h) = Xk[g) Xkih), such 
that the partition factorizes for all g^ and the sum can be 
performed. As a result one obtains the delta function on 
Z 2 , denoted S^^\k), on the vertices of the lattice, which 
require that the representations fee on the edges e meeting 
at the vertex v sum to zero (modulo 2). 

Following this transformation one can define the dual 
theory to the Ising model [67], by introducing new vari¬ 
ables on the vertices of the dual lattice, which depend on 
the {ke} in such a way that all constraints on the original 
vertices are satisfied. Remarkably, the dual variables are 
again Ising spins (elements of Z 2 ), which like the originial 
Ising model only interact with their nearest neighbours. 
The new coupling constant depends on the original /?, 
yet is large if /3 is small and vice versa. Hence, the Ising 
model is also called a self-dual theory, which relates the 
low / high temperature regime of the Ising model to the 
high / low temperature regime of another Ising model 
respectively. 

This particular property of the Ising model will sim¬ 
plify the discussion of the coupling to the gravitational 
background in the following section. 


Xo{g) — 1 Vg G Z 2 , (2.25) 

»<’) = {-! ij'/b-l (2*) 

The Fourier transformed edge weights ujk are then given 
by: 


. _ Jcosh(/3) for fc = 0 
l^sinh(/3) for k = 1 


(2.27) 


C. The coupled system — some basic observations 

In the previous sections IIA and IIB we have intro¬ 
duced two discrete systems, a quantum gravity inspired 
one and a simple matter model. The goal of this work 
is to sensibly couple these models and study their col¬ 
lective dynamics. To be more precise, the gravitational 
part is intended to be the background on which the Ising 
model ‘lives’, providing a dynamical notion of lengths 
which should influence the Ising spin interactions. To do 
so, we work with the following assumptions: 


Eventually, we can expand the partition function into 


Both theories live on the ‘same’ lattice®, in the 
sense that we do not build in a different scaling 
behaviour by hand, e.g. say that the gravitational 


For non-Abelian groups this only applies to class functions, 

which have the property aj(^) = for allE G. Clearly, ^ This can in principle mean that one theory lives on the dual 

all functions on Abelian groups are class functions. lattice, a possibility we do not wish to exclude a priori. 





degrees of freedom should be ‘finer’ than the Ising 
degrees of freedom. This choice is motivated by 
interpreting the lattice as a regulator for both the¬ 
ories, where additionally the lattice itself does not 
carry a notion of length - this notion is supposed 
to be provided by the gravitational theory. 

• The two theories will be coupled by modifying 
the matter part, concretely by implementing a de¬ 
pendence on the gravitational labels into the edge 
weights of the Ising model. One can interpret this 
as a coupling constant /3 depending on the spins of 
the gravitational model, for which we will introduce 
two schemes below. This choice is motivated from 
the coupling of matter to general relativity: While 
gravity is not affected by the matter part in the 
action, the matter part is sensitive to gravity via 
the dynamical metric that enters in the volume 
element, the contraction of spacetime indices and 
covariant derivatives. 

• For the rest of this paper we will consider only reg¬ 
ular (with respect to the combinatorics) square lat¬ 
tices, i.e. we consider 4-valent vertices where every 
pair of vertices share at most one edge. The main 
reason is that coarse graining via tensor network 
renormalization [57, 58] can be straightforwardly 
implemented and will allow us to study the phase 
structure of the model. Moreover the dual lattice is 
also simply a shifted square lattice and we can com¬ 
pare the coupled system to the pure Ising model on 
an equilateral square lattice, e.g. observe changes 
in the phase transition temperature. 

As we have emphasized before, we will interpret the grav¬ 
ity inspired theory as a dynamical background on which 
we place the Ising model. This can in principle be done 
in two ways: 

One possibility is to place the Ising spins on the dual of 
the original lattice, which is demonstrated in fig. 1. In¬ 
deed, a spin foam model coupled to a pure (Yang-Mills) 
lattice gauge theory has been defined on a triangulation 
in [22]. While the 4D Barrett-Crane model [68] ‘lived’ 
on the dual 2-complex, the lattice gauge theory has been 
placed on the triangulation, where the spin foam model 
provided the geometric data, i.e. areas of triangles and 
volumes of 4-simplices. These two theories have been 
coupled by implementing the dependence on the geomet¬ 
ric data into the face weights of the lattice gauge theory: 
via scaling arguments to obtain the proper continuum 
limit, the Wilson action has to scale like the volume of the 
building block, here the 4-simplex, however the plaque- 
tte variables scale like the area squared of the plaquette. 
This mismatch^*^ has to be cured by a proper normali¬ 
sation, which could be interpreted as a locally modified 
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FIG. 1. The coupled model, with Ising spins on the dual 
lattice. The gravitational model is solely placed on the black 
lattice, with SU(2)fe projectors / intertwiners 'Pv{{je}eDv) on 
the vertices v and representations je on the edges. The blue 
dual lattice carries the Ising model, with Ising spins on 
the dual vertices v* and edge weights oje* on the dual edges 
e*. Note that every edge e pierces its dual edge e* only once, 
such that we can interpret the spin je as the length or distance 
between the two Ising spins on e*. There we will implement 
the interaction between the Ising and background degrees of 
freedom, indicated by the blue-black dots. 


coupling constant of the lattice gauge theory, dependent 
on the geometry arising from the spin foam. 

However, the definition of the two coupled theories on 
dual lattices / complexes has a drawback that has been 
pointed out in [24]: Since both theories live on related 
but essentially different lattices, the states and Hilbert 
spaces one defines in the canonical theory, which are es¬ 
sentially product states [8] are not boundary states of the 
spin foam model, such that this model cannot be used 
to compute transition amplitudes between states of the 
canonical theory. To avoid this, one should rather define 
both theories on the same discretisation. In the context 
of our toy model, this version is illustrated in fig. 2. 

In our concrete case, the chosen discretisation of the 
Ising model (with respect to the background) is less im¬ 
portant, since the Ising model is self-dual; essentially we 
explore the same system with inverted temperature^^. 
Nevertheless, these two choices are still worth exploring 
since the choice will impact the coarse graining procedure 
explained in section HI, for which the partition function 
is written as a contraction of (local) tensors assigned to 
the vertices of the lattice. To achieve this in the two dif¬ 
ferent discretisation schemes requires a different choice of 
(equivalent) variables for the Ising spins. In the dual case, 
we stick to the original Ising spins on the dual vertices, 
whereas in the ordinary discretisation it is preferable to 
use the Fourier transformed ones. 


In ordinary 4D lattice gauge theory the underlying lattice is hy- 

percubic and equilateral, such that the area squared of any pla- This is also true for the Ising model coupled to the background 

quette scales as the volume of a 4D hypercube. via a j dependent coupling constant, as introduced below. 
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FIG. 2. The coupled model, with Ising spins on the same 
lattice, after the Fourier transform to irreducible represen¬ 
tations fee has been performed. The gravitational model is 
represented by the black lattice, with SU(2)fc projectors / 
intertwiners 'Pvi{je}eDv) on the vertices v and representa¬ 
tions js on the edges. The blue lattice carries the Ising 
model, with irreducible representation fee and (transformed) 
edge weights on the edges e and a ‘Gauss’ constraint 
■= on the vertices. Again, we 

interpret the 811(2)*, spin js as the distance between two Ising 
spins on the vertices, which enters into the edge weights 
such that the interaction of the two sets of degrees of free¬ 
dom is only via the edges. To underline this, we have drawn 
both lattices distinctly, with blue-black dots indicating the 
interaction on the edges. 


After this extended discussion on the discretisation 
schemes, we have to focus now on the coupling of the 
background and the Ising model, for which a priory many 
choices are possible. As discussed above, we implement 
this coupling by solely modifying the Ising model, i.e. 
by implementing a dependence on the background vari¬ 
ables in the edge weights cOe (or conversely their Fourier 
transforms w*,), concretely by choosing /3 to be a function 
of the underlying geometry, here of the spins je', hence 

P{je)- 

The reader might wonder, why this is supposed to 
be a reasonable choice. Again the reasons stem from 
analogies. First of all, the toy model at hand is sup¬ 
posed to give insights into quantum gravity coupled to 
lattice gauge theory, since 2D spin systems share statis¬ 
tical similarities to 4D lattice gauge theories [28]. If we 
consider again the classical continuum coupling of grav¬ 
ity to Yang-Mills theory, we observe that the metric only 
couples to the curvature of the Yang-Mills connections, 
not directly to the connection. Thus, we rather avoid 
coupling the background degrees of freedom directly to 
the Ising spins, but rather their edge-wise product. In 
such a way we also avoid breaking the global Z 2 symme¬ 
try of the Ising model, such that it can still be Fourier 


transformed and remains self-duah^. Moreover, a back¬ 
ground dependent coupling constant is consistent with 
our intuition of the Ising model. In both discretisation 
schemes, we can interpret the spin je as governing the 
distance between two spins on this (dual) edge. If je is 
increased, the distance between the spins should also in¬ 
crease and result in a weaker interaction, hence a smaller 
effective . Finally, the idea of local coupling constants 
is a notion discussed for lattice gauge theories on irregu¬ 
lar lattices [47] . Here the irregularity of the lattice arises 
from irregular colourings of background spins je, which 
can serve as a guideline for the choice of j3{je)- 

In the sections IIC 1 and IIC 2 we will present two 
different scenarios on how to model the dependence of /3 
on je- 

1. The ‘length’ coupling 

The first approach is a rather simple one, since it only 
takes the distance between two spins into account, thus 
called ‘length coupling’. The coupling looks as follows: 

/3(je) := , (2.32) 

Je+ 2 

where /3o is the standard parameter of the Ising model. 
je + 5 is chosen to avoid an ill behaviour for je = 0. Fur¬ 
thermore it is the edge length appearing in the asymp¬ 
totics of the Ponzano-Regge model [69, 70], a model for 
3D discrete Riemannian quantum gravity. This Ising cou¬ 
pling /3(j) is large for small spins, actually largest for spin 
j = 0. In order to have a direct comparison to the Ising 
model on an equilateral lattice, we rescale /3o, such that 
P'o = Then one obtains the standard Ising coupling 
constant if all j = 0 . 

From the definition of this coupling constant, we ob¬ 
serve already a few basic properties. For /3o ^ 1, i.e. 
high temperature, the Ising model is very insensitive to 
the size of the edge lengths, since all configurations come 
with (almost) the same strength. Hence, for /3o = 0 we 
recover the background without Ising model on it. How¬ 
ever for growing /3o, we notice that the couplings fall off 
with roughly such that we can expect two things: 
First, if je > 0 are excited / allowed, the phase transi¬ 
tion of the Ising model might be pushed towards larger /3o 


Note that we are not implying that a direct coupling breaks this 
symmetry necessarily, but a well-defined definition might not 
be obvious. E.g. a direct product of Yang-Mills and gravita¬ 
tional holonomies might not exist for certain choices of symme¬ 
try groups or just work in a common representation. A direct 
coupling seems to be rather suited if a larger symmetry group is 
used to describe the coupled degrees of freedom. 

Of course, the interpretation of this model is significantly differ¬ 
ent from the standard idea of the Ising model, where the spins 
sit in a rigid lattice (or rather crystal). 
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in comparison to the standard case, since the Ising spins 
effectively see a weaker coupling. Secondly, for /3o S> 1, 
/3(0) » /3(j)Vj ^ 0, such that we expect the Ising model 
to greatly favour j = 0 configurations and thus might 
induce a phase transition in the background from ‘j > 0 
excited’ to ‘only j = 0 excited’. 


8. The ‘area-length’ coupling 

The second option we will investigate is far less ob¬ 
vious, but also interesting. Instead of just assigning 
weaker and weaker coupling constants /? to larger spins, 
we rather want to modify the coupling constant by a 
dimensionless quantity. Therefore we have to compare 
the edge length of the edge to another geometric quan¬ 
tity; the only other local option is the (squareroot of the) 
area the edge is part of. This idea is quite similar to the 
modification done in 4D [22, 23], where one considers the 
ratio of the area squared of a triangle and the 4D vol¬ 
ume of a 4-simplex. In 2D, independent of the chosen 
discretization (or valency of the vertices), every edge is 
shared by two faces, such that we define: 

me}) , (2.33) 

Je +2 

where Ar(/) denotes the area of the face /. 

First of all we note that, if we restrict our discussion to 
square lattices, in which all edges have the same length, 
we recover the standard Ising coupling constant. In the 
case of a dynamical background, we are required to define 
the area of the face / as a function of the edge lengths of 
the face. Instead of constructing a suitable area operator, 
we will go with a poor man’s version: we simply multi¬ 
ply the edge lengths of two adjacent edges, i.e. spanning 
two orthogonal directions, as if we had rectangular an¬ 
gles at each corner of the tetragon. Of course, these will 
give very different answers depending on which pair of 
edge lengths we pick, but this can be improved e.g. by 
averaging over these four choices. 

Clearly, such an area coupling is more non-local than 
the length coupling previously discussed, actually to a 
varying degree. E.g. the area averaging renders the edge 
weight dependent on seven edge labels, which are also 
subject to additional constraints on the bounding ver¬ 
tices of the squares, in particular (generalized) triangle 
inequalities. Therefore, to define this coupling informa¬ 
tion of six vertices of the (background) lattice is required, 
which is obviously at odds with the local coarse graining 
scheme^^. Instead we simplify this interaction by intro¬ 
ducing a more local coupling constant, defined for each 
pair of edge and vertex: 


To disentangle this dependence and write down an amplitude for 
each vertex of the lattice is non-obvious. 


me}) 




i/left 


h + 


bright 



(2.34) 


where jieft and jright denote the edge labels ‘left’ and 
‘right’ of the edge e as seen from the vertex v. In short, we 
simply define the area of the squares locally from the edge 
lengths of the local vertex. Again, if all edge lengths are 
equal, we recover the standard Ising coupling constant. 
Moreover, note that this coupling constant is invariant 
under a global rescaling of all edge lengths. 

In comparison to the length coupling, general state¬ 
ments on the possible behaviour of the composite model 
are not as obvious. Of course, for /3o 1, we recover the 

uncoupled background, since the Ising model essentially 
assigns constant weights to all background configura¬ 
tions. However for larger /3o no immediate consequences 
can be read off, since the coupling on the edge is not 
only determined by the local je, but also its neighbour¬ 
ing ones. E.g. if the neighbouring spins are larger than 
je the coupling constant will be larger than I3q (and thus 
larger than the configuration for all je = 0). But that 
statement is not sufficient to determine whether larger 
spins get suppressed or not, since it only concerns one 
edge: A configuration that increases the coupling con¬ 
stant for one edge generically assigns a lower coupling 
constant to the neighbouring ones. Hence the question 
of possible phases depends on the weights and number of 
certain configurations, yet a general mechanism for sup¬ 
pressing spins larger than j = 0 for /3o 1 does not seem 

to exist. 

To make this statement more precise, let us consider 
a configuration, which is allowed by the coupling rules 
and deviates from the regular square case. Eor the time 
being, we ignore the amplitude of the background and 
just consider the influence of the Ising edge weights. A 
typical configuration is e.g. twice spin j = 0 and twice 
spin j = j' > 0 on the four edges - all possible permu¬ 
tations of these spins are allowed by the coupling rules 
of SU(2)fc. If we then consider the product of all four 
edge weights, e.g. for alternating spins and all ke = 0 
normalized by the equilateral case, we obtain: 


ni=l ^oijiT ji+l ^ ji—l) 

(wo(0,0,0))4 


cosh^ 



cosh'‘(/3o) 



(2.35) 


since the spins j' are ‘surrounded’ by spins j = 0 and vice 
versa. Thus we have a hyperbolic cosine squared with 
an effectively larger coupling and one with an effectively 
smaller one. Similarly this also occurs for more general 
cases with non-equal spins {ji} on the four edges of the 
vertex. 

Remarkably, it turns out that the product of four edge 
weights (for the same kg configuration) for configurations 
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deviating from the all spins j equal case is larger in gen¬ 
eral than for equilateral configurations. Thus, the Ising 
spins effectively see a larger coupling if the lattice spac- 
ings are not all equal. The explanation is the exponen¬ 
tial growth of the hyperbolic cosine (or sine for kg = 1), 
where the shift to larger /3o dominates over the shift to 
smaller (still positive) Pq. Therefore, this weight is par¬ 
ticularly large if the four spins j around a vertex are 
very different: The maximal effective coupling occurs if 
the spins on the edges of the vertex alternate between 
smallest {j = 0) and largest spin of the quantum group 
U = jmax) “ the most ‘non-equilateral’ configuration, see 
(2.35) for j' = jmax- If the deviation from the equilateral 
case is very small, e.g. only spins j and j -I- 1 are excited 
with j ^ 1, then the effective coupling quickly converges 
to the equilateral one. 

Thus we can expect a very different behaviour from 
the ‘area’-coupling than the ‘length’-coupling: A mech¬ 
anism to suppress configurations with larger spins j, ap¬ 
pears to be absent in the ‘area’-coupling. Preferred con¬ 
figurations are rather characterized non-locally by their 
(ir)regularity, where the lowest weight is assigned to the 
equilateral configurations (independent of the spin j) and 
the highest weights are assigned to the most irregular 
configurations alternating between smallest and largest 
spins. One implication is that configurations with spins 
j ^ 0 will occur that are equipped with a larger weight 
than the configuration with all j = 0. Thus one can ex¬ 
pect a shift of the Ising phase transition towards smaller 
Pq if more spins j are allowed, since the coupling is ef¬ 
fectively larger. A second implication concerns the back¬ 
ground, which might be affected by the preference of ir¬ 
regular configurations^®. 

In order to answer these questions and properly exam¬ 
ine the differences between the two couplings of the Ising 
model to the background, it is imperative to study the 
dynamics of the system beyond just a single vertex. To do 
so, we will introduce a coarse graining algorithm, known 
as tensor network renormalization, in the following sec¬ 
tion, which will allow us to extract the phase structure 
of the coupled system. 

III. COARSE GRAINING ALGORITHM - 
TENSOR NETWORK RENORMALIZATION 

In section II we have thoroughly motivated and intro¬ 
duced a system of Ising spins coupled to a dynamical 
background as a toy model for lattice gauge theory cou¬ 
pled to discrete quantum gravity, where both theories are 
regularized by the same (or the dual) lattice. The goal 
of the remainder of this work is to study the composite 


For /3o —> oo, this weight for irregular configurations diverges 
with respect to all j = 0. Since the background always couples 
to the trivial representation and is normalized with respect to it, 
the background probably will not stay unaffected. 


dynamics of this novel system. By this we particularly 
mean the composite phases of the system, i.e. the regions 
in the parameter space of the theory, in which the system 
shows the same qualitative behaviour. An example would 
be the Ising model in the (dis) ordered phase on a specific 
background, e.g. all spins j < J are allowed. Moreover, 
we intend to investigate the quantitative changes each 
of the models experiences, e.g. whether the phase tran¬ 
sition temperature of the Ising model shifts due to the 
deviations of the background away from equilateral lat¬ 
tices and whether the background is affected by the Ising 
model. 

In this context, the discretisation plays a crucial role. 
Unless we discuss topological theories, the results of the 
theory will depend on the chosen discretisation. In the 
standard context, this discretisation inherits a length 
scale from the fixed background and one describes how 
the dynamics of the system changes as one considers a 
different scale, which is the standard idea of renormaliza¬ 
tion. Since a background scale is absent in our situation, 
we have to relate theories defined on different discreti¬ 
sations instead, e.g. the system on a finer discretisation 
to a system on a coarser discretisation, via a transforma¬ 
tion, here coarse graining, which should not change the 
partition function and expectation values of observables. 
We can understand this transformation as a map from 
finer to coarser degrees of freedom. 

A coarse graining algorithm well suited for studying 
systems without a direct reference to a length scale is ten¬ 
sor network renormalization [57, 58], originally developed 
in condensed matter physics. Note that Monte-Carlo 
simulations are not applicable, since the intertwiner mod¬ 
els [59] are inherently complex; the same also holds for 
spin foam models. The basic idea of this algorithm is 
to rewrite the partition function as a contraction of mul¬ 
tidimensional arrays, the tensors, which encode all the 
dynamical information of the system^®. Frequently the 
tensors have the same number of indices as the original 
lattice and are pictorially represented by a vertex with as 
many legs as indices. An index contracted with another 
tensor is represented by connecting the two legs of the 
tensors. Hence the partition function is represented by 
many tensors connected to one another according to the 
combinatorics of the lattice - a network. 

The range of the index is frequently referred to as the 
bond dimension y. In our case, as we deal with finite 
groups and quantum groups that come with a natural 
cut-off on the representation labels, the bond dimension 
of the (initial) tensor is finite. However a straightfor¬ 
ward definition for systems with (compact) Lie groups. 


For lattice gauge theories, tensor network representations exist 
as well, see e.g. [51, 71], but are very costly due to storing re¬ 
dundant (gauge) information. In [50] an algorithm for 3D lattice 
gauge theory has been implemented and tested, in which not all 
information is stored in the tensor network. The data encod¬ 
ing gauge symmetries are rather used to ‘decorate’ the tensor 
network, with the gauge conditions solved. 
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FIG. 3. The initial amplitudes for the dual and the standard 
discretisation. The partition function is chopped into local 
amplitudes, tensors, along the edges. Since the same vari¬ 
ables are shared between the separated tensors, only the edge 
weights, tUe* and respectively, are split. Then the ampli¬ 
tudes are combined into tensors, a ‘decorated’ tensor for the 
dual theory, and a ‘standard’ tensor for the normal one. 


e.g. U(l) or SU(2), as their underlying symmetry group 
is not obvious, unless one can introduce a approxima¬ 
tive cut-off as in the strong coupling expansion or use a 
different scheme to write the system with a finite bond 
dimension [72, 73] or incorporate Monte-Carlo methods 
for the contraction of indices [74]. 

In our concrete example, we chop the partition func¬ 
tion into local amplitudes, which we can associate to the 
vertices v of the lattice. Schematically, we show this in 
fig. 3. These local amplitudes are summarized as a tensor 
T on the vertex v of the lattice: 


'Pv{{je}eDv) ^i9s{e*)9t{e*) i {je}) j 


e*dual to eZ)v 


(3.1) 


= VA{Je}e^v) 6^^^ ( • 

\eZ)v / eZ)v 


(3.2) 


Note that there is a slight difference in the definition 
of the tensor if we have chosen either the lattice or its 
dual for the Ising spins. In the former case, the standard 
construction goes through, since the variables solely live 
on the edges e and their mutual dependence is encoded 
on the vertices v. In the latter case, the dual, this only 
applies to the background, however the Ising spins live on 


the dual vertices v* with edge weights on the dual e*. In 
fact, splitting of the edge weights is straightforward, but 
one has to keep track of the Ising spins, a ‘decoration’ to 
the tensor network. How to coarse grain such ‘decorated’ 
tensor networks has been developed in [50] and turns 
out to be as straightforward as the standard algorithm 
[57, 58]. Due to the self-duality of the Ising model, we 
nevertheless drop this option and focus on the model, in 
which both systems live on the same lattice, see also fig. 
2. To simplify the notation, we will from now on drop 
the superscript. 

To obtain the whole partition function again, the ten¬ 
sors are connected, i.e. the shared indices are summed 
over, and the partition function is written as the tensor 
trace of the tensor network: 

Z = Ttr n(r T...) . (3.3) 


So far, we have obviously only rewritten the problem into 
a different form, which is more local than the original 
one^^. The idea of the algorithm is then to locally trans¬ 
form / coarse grain the tensor network, such that origi¬ 
nal partition function is approximated by a coarser, still 
local tensor network. In short, the tensors themselves, 
and as such the dynamical ingredients of the system, get 
renormalized. 

There exist a plethora of different tensor network al¬ 
gorithms [57, 58, 75-77], which more or less only differ 
in the scheme of how several fine tensors are transformed 
into a coarse one. Such a transformation can be straight¬ 
forward, e.g. on a square lattice, one can simply combine 
four tensors on the corners of a square into one by con¬ 
tracting the inner legs. One straightforwardly arrives at 
a new tensor, however with a squared index size, pic- 
torially represented by double edges. Continuing with 
this endeavour, also with respect to inevitable numerical 
simulations, is fruitless and requires a truncation for two 
reasons: a practical one, since no computer possesses infi¬ 
nite memory and infinite computational time to contract 
the indices, and an interpretative one, in order to com¬ 
pare the coarse to the fine tensor. Hence one is required 
to define new coarse degrees of freedom arising from the 
fine ones, such that one preserves an interpretation and 
can simultaneously truncate the number of degrees of 
freedom with a good control on the error being made. 

At this stage, we will not repeat the whole introduc¬ 
tion and derivation of tensor network renormalization 
and instead briefly discuss the triangular algorithm, orig¬ 
inally introduced in [50] . The derivation of the formulae 
is straightforward from the related 4-valent algorithm 
[54, 57, 58], which we explain in more detail in appendix 
B. 
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Here local means that each tensor is only connected to a nearest 
neighbour tensor. 
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A. The triangular algorithm 

As the name suggests, the triangular algorithm shifts 
the perspective from 4-valent tensors to 3-valent ones: 
Instead of starting with one tensor T, we start with four 
tensors Si. These tensors can be obtained e.g. from 
the first step of the 4-valent algorithm, in which one 
splits the 4-valent tensor across both diagonals into two 
pairs of 3-valent ones, see appendix B for a more detailed 
explanation. 

The triangular algorithm itself is shown in fig. 4: To 
compute a new 3-valent tensor, two of the Si are com¬ 
bined pairwise, where the shared edge gets contracted, 
resulting in a new 4-valent tensor. To arrive again at a 
3-valent vertex, the two ‘parallel’ edges in fig. 4 have 
to be combined into one, since they represent the ‘finest’ 
degrees of freedom of the system; for this we have to 
introduce an embedding map, which is computed via a 
singular value decomposition (SVD). 

Let us discuss this for a more concrete example: Sup¬ 
pose we want to construct the S[ in fig. 4 from S 2 and 
S 4 . We define the intermediate tensor Si as: 

{Si)ab-jl = ^(S'2)o;fci iS4)b-jk ■ (3.4) 

k 

We have to define an embedding map for the indices {jl) 
into a new effective index c'. Therefore we perform a 
SVD on the following matrix: 

CSi\abym=ll^ia>.)4^^Vui), ■ ( 3 . 5 ) 

i 

U(ab). i and V(j 7 ) i are the singular vectors, ^ the singular 
values of the matrix iSi)(a.b)-,(ji)^ where the indices (ad) 
and (jl) simply denote product indices. U and V are uni¬ 
tary matrices, i.e. the singular vectors are orthonormal, 
while the singular values are non-negative and ordered 
in size, Ai > A 2 > ... > Aat > 0 . 

U and V can be interpreted as mapping the indices (ab) 
and (jl) to i respectively. The (relative) value of Ai, e.g. 
with respect to the largest singular value, determines how 
important the index i is in reconstructing Si. Thus we 
use V both as an embedding map, mapping the indices 
(jl) into a new effective index i, and as a truncation, 
cutting off the index range of the new edge. To do so, we 
have to multiply it from the right to the tensor Si and 
sum over the indices (jl). Then we obtain the new S[ as: 

{S'i)c’-,ab = U(ab),c'^c' , (3.6) 

where we have used the fact that V is a unitary ma¬ 
trix. However, in order to leave the partition function 
unchanged, it is necessary to contract the new S'l on the 
other side with obtained from Si, as it is also im¬ 
plied in fig. 4. Hence the formula for S 2 contains the 
explicit contraction^®. Also note that this algorithm is 


In fact, it is necessary to compute the embedding maps from 



J I 



S[ 


c' 

FIG. 4. The triangular algorithm. Instead of one 4-valent 
tensor, forming a rectangnlar tensor network, we work with 
fonr 3-valent tensors Si. We glue them pairwise into inter¬ 
mediate 4-valent tensors, compnte an embedding map from 
this tensor, here called U and V, which maps two fine edges 
into one effective coarse edge withont affecting the partition 
function since UU'^ = = 1. Thus we obtain four new 

3-valent tensors S',' and iterate the procedure. 

not as ‘symmetric’ as the 4-valent one: Every 3-valent 
tensor has a special edge, namely the one obtained after 
applying the latest embedding map. 

Of course, similar to the 4-valent algorithms employed 
in [51-54], it is beneficial to explicitly preserve the un¬ 
derlying symmetries of the tensors under coarse graining 


both Si and S 2 and check, which one of them leads to a smaller 
overall error and thus a better approximation, similar to the 
higher order SVD used in [76]. 
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because of saving of computational resources and keeping 
an interpretation of the new coarse degrees of freedom in 
terms of the original variables. Thanks to representation 
theory, we know that some couplings are forbidden, e.g. 
for the Ising model if the Gauss constraints are not sat¬ 
isfied. To exploit this as much as possible, one writes the 
matrices that are to be split in a block diagonal form; 
we call this the recoupling basis. In the non-Abelian or 
quantum group case, this amounts to a precontraction 
of magnetic indices, as the dependence on those is not 
changed under coarse graining. Each block is then la¬ 
belled by a set of representation labels, e.g. (j, k) for our 


toy model, called the ‘intertwiner channel’. Then instead 
of performing the SVD on the entire matrix at once, one 
performs one SVD per block matrix, such that one en¬ 
dows the block label onto the new effective edge. 

For the triangular algorithm, this symmetry protecting 
algorithm can be straightforwardly derived from its 4- 
valent counterpart. Therefore, we here only present the 
formulae, the actual derivation of the 4-valent version is 
put in appendix B. 

Let us start with the block diagonal form of the 3- 
valent tensor, which is similar to equations (B 8 ) to (Bll) 
of the 4-valent algorithm: 


32 3i 


h) Y ’ 

(3.7) 

45 


kb) j' , 

(3.8) 


J3 J4 

h 

(^ 3 F ^^({ n { 2 . 3 }, WI 2 . 3 }) =( 4 )(^^’'=^)({ n { 2 . 3 }) <^^^^*2 + ^3 - k ,) , 

J3 

jl 


(3.9) 


(3.10) 


J4 


where we have summarized the representation labels 
ke as {le}- The diagrams encode the dependence on the 
magnetic indices, which remains unchanged under coarse 
graining due to preserving the symmetries. 

Analogue to the 4-valent algorithm, the next step is 
the computation of the intermediate 4-valent tensor in 
its block-diagonal form to which the SVD will be applied. 


Remarkably, this expression can almost immediately be 
read off equations (B12) and (B13), together with the 
fact that identity (B14) splits into two 6j symbols. This 
allows us to essentially cut equations (B12) and (B13) in 

half to obtain the equations for and S 2 - 


J-Os.fes) 


Jr(is,*5) 


(A. 4, /.) = V V ysyi- 


Jc Ja 
jl 32 


b {m} 


(A. A; A. A) = V V 72-3- 


Jc Ja 
J4 ja 


d {m} 


3b 

3b 


S^^\ki + k2 - k^) 


ja 

jd 


d^^Hks + k4 - k^) 


(3.11) 


(3.12) 


Analogously, one obtains the equations for S 3 and S '4 from the equation of T 2 , but essentially it is merely a suitable 
permutation of labels in the equation. Eventually, we can compute the embedding maps by performing a SVD on Si, 
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e.g. on 


V /(/l,/2);(/c./a) 


E ttUuM) \{ 35 M) 




{Ia,Ia.),^l 


il 


(3.13) 


where ii denotes the multiplicity of the channel (jsjfcs) in this particular example. If this is the proper embedding 
map, i.e. the one with least error approximating both for Si and S 2 , then we eventually obtain for the new 


[S'l)^- 




(ST\, ,_ rrOs.fcs) \05,fc5) 


(3.14) 


whereas we find for S' 2 : 


(^'2)(' 


jbMM 


'({n{3.4}) = E^2 


_j^)(yt)0a.fc5) 


(/c./a).n 


(3.15) 


The derivation of S'z and ^'4 works analogously. 

Fortunately, in the models under discussion here, all 
the Si turn out to be identical, also under coarse grain¬ 
ing because the same recurrence relation holds for all of 
them^*^. As a result, it will be sufficient from now on 
to just work with one 3-valent tensor instead of four, 
which simplifies equation (3.11), but more importantly 
always allows us to use (3.14). Thus we greatly reduce the 
amount of data necessary to be stored during the numer¬ 
ical simulations and additionally avoid the explicit con¬ 
traction of (roughly) 2 x 2 x jmax indices. Here this might 
appear to be insignificant, but as soon as more compli¬ 
cated systems are examined, for example SU(2)fc x SU(2)fc 
spin nets, which roughly require quartic amount of re¬ 
sources in comparison to the background discussed here, 
these gains are crucial. 


1. Truncation scheme and identification of phases 

For our study of the presented toy model, we will use 
a very simple truncation scheme for the singular values. 
Instead of comparing the singular values obtained from 
all blocks and keeping the y largest of these, we take only 
one singular value for each block (j, k). This certainly is 
a drastic simplification, but it seems sufficient to capture 
many important features of the model: Taking one singu¬ 
lar value for each block is close to the initial definition of 
the models, in particular we can capture the fixed points 
of both the background and the Ising model and thus can 


We assume that the embedding maps obtained from Si is the 
one giving the smallest error, therefore we also suppress the ad¬ 
ditional superscript (1). 

This is due to the choice of splitting of the signs in (B13) sym¬ 
metrically in (3.11) and (3.12). In principle, a different splitting 
is possible, e.g. assigning the sign just to one of them. However, 
we have found that our choice is numerically more stable and 
consistent with results of the 4—valent algorithm. 


qualitatively study the phase diagram of the model. Cer¬ 
tainly, we are aware that this approximation breaks down 
in certain regimes of the model, e.g. close to the phase 
transition of the Ising model, where more singular values 
have to be taken into account, which results in shifts of 
the position of the phase transitions. However, numeri¬ 
cal precision is a rather secondary concern, since we are 
dealing with a toy model tailored to demonstrate the po¬ 
tentially rich dynamics of matter coupled to spin foams. 
Moreover, as many other tensor network renormalization 
methods, our scheme requires an infinite bond dimension 
to study (second order) phase transitions. A recently 
developed algorithm [77] allows the study of phase tran¬ 
sitions at relatively low bond dimension - we leave an 
adaptation of our current algorithm for future research. 

Fortunately, this truncation scheme allows for a 
straightforward identification of the phases of the model 
via the singular values in the intertwiner channels (j, fc), 
j for the background, k for the Ising model. These 
can then be summarized in a (jmax -|- I) x 2 matrix. 
E.g. the topological background given by J (see again 
section HA) is given if all intertwiner channels j < J 
have one non-vanishing singular value equal to one and 
all channels j > J have only vanishing singular values. 
Hence the singular values directly indicate which spins 
j are allowed. Similarly for the Ising model, the dis¬ 
ordered phase is given if only the channel k = 0 has 
one non-vanishing singular value equal to one, the or¬ 
dered phase if both k = 0 and k = 1 have one non¬ 
vanishing singular value equal to one. Let us conclude 
with an example for the combined model: If only the 
channels {j,k) = ( 0 , 0 ), ( 1 , 0 ), ( 2 , 0 ), ( 0 , 1 ), ( 1 , 1 ), ( 2 , 1 ) 
possess one singular value equal to one while all others 
vanish, this indicates the background model J = 2 and 
the ordered Ising phase. 

In the sections IV and V we will discuss the results 
obtained from coarse graining the Ising model coupled 
to the dynamical background. In the next section IV 
we first focus on the influence of the background on the 
Ising model by choosing a topological fixed point for the 
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background. 

IV. THE ISING MODEL COUPLED TO A 
TOPOLOGICAL BACKGROUND 

In the remainder of this paper, we will present the 
results of applying the coarse graining algorithm intro¬ 
duced in section III to the model of Ising spins coupled 
to a dynamical background discussed in section IIC. As 
we have explained in section II A, there exist many dif¬ 
ferent choices for the background, encoded in the factors 
a({j}) modifying the Clebsch-Gordan coefficients of the 
quantum group. Recall that these factors are chosen such 
that the local amplitudes are invariant under all 2D Pach- 
ner moves, such that they describe topological systems. 
Additionally, one can also consider linear combinations of 
the factors a({j}) to go beyond topological backgrounds. 

To order the results pedagogically, we start by dis¬ 


cussing the cases in which the background is given by 
a topological theory, that is only one a({j}) is chosen 
instead of a superposition, which will be addressed in 
section V. Still, this ansatz allows for an interesting inter¬ 
pretation. Since the background by itself is topological, 
it is also a hxed point of the coarse graining algorithm: 
usually for small deviations around this fixed point, the 
system flows back to this fixed point. Here, the Ising 
model at /lo > 0 serves as a deviation, such that the 
composite system is not a fixed point and actually flows 
under renormalization. Hence, it is foremost interesting 
to discuss how this new background geometry affects the 
Ising model, in particular the change of its phase transi¬ 
tion temperature, e.g. if larger edge lengths are allowed. 
For larger /3o we can also expect changes in the back¬ 
ground. 

Let us introduce the initial tensors for the start of the 
coarse graining procedure. The 4-valent tensor T is of 
the general form (see also appendix B, e.g. (B6)): 


J, {fcj) = ^ ({j J, (fci + 

is 


- fcs) 



J3 J4 


(4.1) 


The block diagonal part Ti is given as follows: 




a(jl,j2, j5)o(j3, js) 


\/wfci(ji; {je})Wfc2 O2; {je})Wfc3 (ja; {je})i^kAk\ {je}) 


(4.2) 


Vs 


where the argument of Cj indicates that the Ising edge weights can also depend on neighbouring spins jei e-g- in the 
‘area coupling’. To obtain from (4.2) the amplitudes for the triangular algorithm is straightforward. In case the Cjk^ is 
only of function of je, (4.2) factorizes and can be trivially split; if it is a more non-local function, one has to apply a 
SVD on (4.2) completely analogous to the 4-valent algorithm described in appendix B. The 3-valent tensors obtained 
from that give the initial amplitudes. 


In the following, we will summarize the results first 
for the ‘length coupling’, followed by the ‘area coupling’. 
We discuss the choices of topological hxed points for the 
background presented in section HA for different levels 
k of the quantum group SU(2)fc. 

A. ‘Length coupling’ — Results 

Let us briefly recap the length coupling introduced in 
section HCI. A dependence of the Ising edge weight 
is implemented by modifying the coupling constant 
~ that is the larger the distance between two spins, 
the weaker their interaction. 

From this coupling mechanism, we already conjectured 
that for /3o > 1 spins j > 0 will get suppressed by the 
Ising edge weights and there should exist a phase for 
which the Ising model is in its ordered phase (low tem¬ 


perature) on a background with all j = 0. Conversely, 
in the limit /3o —>■ 0, i.e. high temperature, we reobtain 
the pure background, such that we conjecture that there 
exists a phase with the disordered Ising model on the 
initially chosen topological background. 

All spins J < jrnax allowed: The largest class of 
topological fixed points we are considering for the back¬ 
ground are those which allow all spins j < J, where J € N 
and J < jmax = §• That is for larger level fc, we can con¬ 
sider more models. The results are summarized in tables 
I and H, since we hnd two phase transitions for almost 
all backgrounds: 

First of all, we can verify both previously made con¬ 
jectures: For /3o ^ 1 there exists an extended phase, in 
which the geometry is given by the chosen topological 
background (all j < J allowed) and the Ising model is 
found in its disordered / high temperature phase. For 
/I 3> 1 we find another extended phase, in which the 
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Transition in /^“interval 

fc = 4 

fc = 5 

k = 6 

k = 7 

fc = 8 

J = 0 

[0.388,0.389] 

[0.388,0.389] 

[0.388,0.389] 

[0.388,0.389] 

[0.388,0.389] 

J= 1 

[0.694,0.695] 

[0.716,0.717] 

[0.729,0.73] 

[0.738,0.739] 

[0.744,0.745] 

J = 2 

[0.821,0.822] 

[0.913,0.914] 

[0.969,0.97] 

[1.006,1.007] 

[1.032,1.033] 

J = 3 

n.a. 

n.a. 

[1.061,1.062] 

[1.158,1.159] 

[1.224,1.225] 

J = 4 

n.a. 

n.a. 

n.a. 

n.a. 

[1.295,1.296] 


TABLE I. Interval for parameter /3 marking the phase transition of the Ising model from the disordered to the ordered phase 
in the ‘length’ coupling. The background model, labelled by J, does not change. 


Transition in /^“interval 

fc = 4 

k = 5 

fc = 6 

fc = 7 

fc = 8 

J= 1 

[1.658,1.659] 

[1.693,1.694] 

[1.715,1.716] 

[1.731,1.732] 

[1.742,1.743] 

J = 2 

[1.69,1.691] 

[1.761,1.762] 

[1.81,1.811] 

[1.845,1.846] 

[1.87,1.871] 

J = 3 

n.a. 

n.a. 

[1.838,1.839] 

[1.898,1.899] 

[1.943,1.944] 

J = 4 

n.a. 

n.a. 

n.a. 

n.a. 

[1.965,1.966] 


TABLE 11. Interval for parameter /? marking the phase transition of the background from < J to je = 0 in the ‘length’ 
coupling. The Ising model stays in the ordered phase. 


background is restricted to j = 0 and the Ising model is 
in its ordered phase. However, these two phases are not 
directly separated by a phase transition (of geometric and 
Ising degrees of freedom). In between there exists a phase 
parametrized by the ordered Ising model on the initially 
chosen topological background. Hence we observe two 
phase transitions, the first (for rising ,0o) is a transition 
purely of the Ising model on the topological background, 
whereas the second is a transition purely of the geome¬ 
try, triggered by the particular coupling (via edge lengths 
alone) of the Ising model to the background. 

Besides the appearance of three different phases of the 
composite system, it is interesting to study the influ¬ 
ence of the background system on the position of the 
phase transitions. The case J = 0, essentially the Ising 
model on a regular square lattice, is our reference point. 
Due to the approximation discussed in section HI A I of 
the used algorithm, we find the critical temperature at 
0.388 < Po < 0.389 instead of Pent ~ 0.4406.... This 
also implies that the other phase transition locations are 
subject to change if the accuracy is increased. Hence 
we consider a very precise determination of the phase 
transition temperatures (in this approximation) to be a 
fruitless endeavour and rather interpret the findings more 
qualitatively. 

Let us consider the phase transition of the Ising model 
(see table I) first: We can study this data from two per¬ 
spectives. Either we consider a particular background J 
and study the change of the phase transition if we in¬ 
crease the level k of the quantum group (a row) or we 
examine how the temperature changes for the same level 
k, while increasing the maximum allowed spin J (a col¬ 
umn). 

In the former case, we observe that for a given J, the 
phase transition parameter increases for growing fc, yet 


this increase is not uniform for increasing k. The in¬ 
crease is generically largest between J = jmax and the 
next largest level k, and appears to decrease under in¬ 
creasing k further, such that it might eventually con¬ 
verge. The interpretation of this behaviour is that, for 
the same underlying geometry, the Ising model effectively 
has a weaker coupling /?, such that its phase transition 
occurs for larger Pq, even though the edge weights are not 
changed in different quantum groups. Hence the origin 
of this behaviour must lie in the background geometry: 
Take e.g. the model J = 2 for fc = 4, i.e. J = jmax = 2. 
From the coupling rules / Clebsch-Gordan coefficients, 
one can determine that there exist only three possible 
ways how two spins ji, j 2 can couple to a third spin j = 2, 
namely (0,2), (1,1) and (2,0) - j = 2 here is of quan¬ 
tum dimension one. If we go to fc = 5, there are then 
five ways, namely also (1, 2) and (2,1) in addition to the 
previous ones. As a result, there exist now more configu¬ 
rations in the partition function with larger edge lengths 
and hence an effectively lower temperature P of the Ising 
model, pushing the phase transition towards higher Pq. 
If the quantum group level k is increased further, even 
more couplings are allowed, however these involve spins 
j > J, which are disallowed by the geometry (or only ex¬ 
cited briefly under coarse graining), such that the phase 
transition temperature is less affected. 

On the other hand, if we keep k fixed and increase 
J, we also observe an increase in the phase transition 
parameter Pq, where the gap between J and J -I- 1 de¬ 
creases as one approaches J = jmax of the quantum 
group, the jump from J = 0 to J = 1 being the most 
significant. This is a very much expected result consid¬ 
ering that P{j) falls off with 1/j in the edge weights of 
the Ising model and the fact that increasing J allows not 
only larger spins j, but also more conhgurations of the 
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geometry with non-vanishing j on it. The latter effect 
can be seen e.g. for k = 8 : even though only for J = 4 
the largest edge lengths get excited, the increase in /3crit. 
is rather small. Since j = 4 is of quantum dimension one, 
only few new couplings become allowed with respect to 
the J = 3 model. Thus we see that the weakening of the 
Ising coupling is not solely determined by the (larger) 
edge lengths but also the amount of configurations, deter¬ 
mined by the background, that allow these edge lengths 
to occur. 

For the second phase transition (see table II) , that is a 
geometric phase transition from a topological fixed point 
to a geometry with only j = 0 allowed (with the Ising 
model in its ordered phase in both cases), we perform a 
similar analysis: If we keep J of the background fixed and 
increase the level k of the quantum group, we again ob¬ 
serve an increase in the transition parameter /3o. The 
reason is similar to the previous case: If the level of 
the quantum group is increased, more couplings (of non¬ 
vanishing spins j) become allowed increasing the number 
of configurations with non-vanishing spins. Thus Pq has 
to be increased further for larger k to sufficiently suppress 
spins j > 0 such that the system flows to a geometry with 
all J = 0. A similar argument is valid in case of keep¬ 
ing k fixed and increasing J, where we also increase the 
number of configurations with non-vanishing j. 

This second phase transition is of course a very pe¬ 
culiar example of how the matter system influences the 
background. If the coupling Pq of Ising spins becomes too 
strong, it forces the background to only allow its shortest 
possible distance between two vertices. This is a feature 
one would rather like to avoid for 4D gravity plus mat¬ 
ter, such that one can see this geometric transition as an 
indication that Po should be restricted to smaller values 
or that the coupling P{j) should be modified to favour 
spins J > 0. 

Only j = 0 and j = jmax allowed: The next class 
of models we tested, is the background in which only 
representations with quantum dimension dj = 1 are al¬ 
lowed, that is only j = 0 and j = jmax for even levels 
k of the quantum group. Again, we observe two phase 
transitions; their locations for various k is given in table 
III. 

The phase transitions are identical in type as for the 
previous background. On the one hand we find one tran¬ 
sition at lower /?o, which is a pure transition of the Ising 
model on the background, from the disordered to the or¬ 
dered phase. On the other hand we observe again the 
background transition induced by the length coupling, 
which forces the background to take its smallest edge 
length possible at larger Pq. 

Concerning the Ising transition, the actual position of 
the phase transition is quite different from the j < J 
backgrounds. In general, the phase transition can be 
found at lower Pq (unless only j = 0 is allowed) and in¬ 
creases only slowly as k is increased. This is due to the 
peculiarity of this model: for even k, both j = 0 and 
j = jmax have quantum dimension dj = 1; the only 


allowed non-trivial coupling rule is {jmax, jmax) —> 0. 
Thus, even though the Ising model is added in a non¬ 
trivial way, the background cannot flow away from this 
general configuration, i.e. no other spins j can get ex¬ 
cited during the renormalization procedure. Therefore, 
there are much fewer configurations of the background 
with non-vanishing spins, even in contrast to the low¬ 
est j < J models. Moreover, this remains unchanged as 
k is increased, yet the interaction between neighbouring 
Ising spins on an edge carrying jmax is weaker for larger 
fc, pushing the transition to larger Pq. 

For the geometric transition to a background where 
only j = 0 is allowed, we find a different situation with 
respect to the j < J models. Instead of increasing Pq for 
growing J or fc, we observe a decrease. Again, this is due 
to this peculiar model, in which no other spins except 
j = 0 and j = jmax are (and stay) excited. Since jmax 
grows as k gets increased, the edge weights for jmax fall 
off faster with respect to j = 0 for larger k, resulting in 
a transition to ‘just j = 0’ at lower Pq. 

Only j = 0, j = and j = jmax allowed: The 
last background model to investigate together with the 
length coupling is the background with j = 0, j = 
and j = jmax excited. Since such a background can only 
exist if the level k of the quantum group is a multiple of 
four (this is necessary, but not sufficient, see also section 
II A), we study it for only three cases. The results are 
summarized in table IV: 

As we have also observed for the other backgrounds, 
there are two phase transitions in the length coupling: 
the first one of the Ising model from ‘disordered’ to ‘or¬ 
dered’, the second one in the geometry from the back¬ 
ground to the background with only j = 0. Across 
the different quantum groups, we observe that the Ising 
transition occurs at larger Pq, since the spins j = 
and j = jmax increase for larger k, while the couplings 
among the spins j (and thus the background configura¬ 
tions) remain unchanged. This results again in an effec¬ 
tively weaker interaction between the Ising spins as k is 
increased. For the second transition, we see a decrease 
in the critical Pq from A: = 4 —>■ fc = 8, similar to the 
(0, jmax) background. 

More interesting is the direct comparison across the 
backgrounds (for the same quantum group) because we 
can directly study the changes that more excited spins 
have on the positions of phase transitions. Note that we 
only compare the results of fc = 4 and k = 8 . 

Let us start with the pure Ising transition, while com¬ 
paring the models (0,jmax) and (0,^,jmax)- For the 
latter, we generically observe the transition at larger Pq, 
which is clear, since more configurations with weaker 
Ising spin interactions (with respect to j = 0) are al¬ 
lowed. The relation to the j < J models is less obvious: 
For k = A and J = 1, jmax is not excited, such that 
edges with the weakest Ising spin interactions are forbid¬ 
den, yet the transition still occurs for a larger Pq with 
respect to the (0, jmax) model. This is the case because 
the J = 1 model permits more configurations with non- 
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Transition in bAnterval 

k = 4 

k = 6 

k = 8 

Ising ‘disordered’ —^ ‘ordered’ 

[0.639,0.64] 

[0.666,0.667] 

[0.681,0.682] 

U = 0; j = jmax) j = 0 

[1.25,1.251] 

[1.175,1.176] 

[1.139,1.14] 


TABLE III. Results for the background with only j = 0 and j = jmax allowed for even k together with the Ising model in the 
‘length’ coupling. The first row is the first transition we find, a pure transition in the Ising model form ‘disordered’ to ‘ordered’ 
as /3o is increased. The second transition is geometric, where for large /3o just j = 0 is allowed. 


Transition in /i-interval 

II 

00 

II 

Ising ‘disordered’ —>■ ‘ordered’ 

[0.821,0.822] 

[1.093,1.094] 

(j = 0; j = j = jmax) j = 0 

[1.69,1.691] 

[1.576,1.577] 


TABLE IV. Results for the background with only j = 0, j and j = jmax allowed for fc = 4/, 1 € N, together with the 

Ising model in the ‘length’ coupling. The first row is the first transition we find, a pure transition in the Ising model from 
‘disordered’ to ‘ordered’ as j3o is increased. The second transition is geometric, where for large /3o just j = 0 is allowed. 


vanishing spins, giving less weight to the ‘only j = 0’ 
case. The J = 2 model is (up to signs) identical to the 
(0, jmax) niodel, both show consistent results. Since 
they allow all configurations of the background allowed 
by the quantum group, they show the largest phase tran¬ 
sition parameters for k = A. For k = 8 we observe a sim¬ 
ilar behaviour across the different backgrounds: despite 
the small weight assigned to edges carrying the maxi¬ 
mal spin jmax, the transition occurs first for the model 
(0,jmax) followed by the j < J models for J = 1 and 
J = 2. Remarkably, the latter is already very close to the 
transition in the (0, jmax) model, /3o gets increased 
further for J > 2. Again we conclude that the amount 
of background configurations with non-vanishing spins 
j quickly outweighs the weakening of Ising spin interac¬ 
tions by large spins j. 

The observation is even more striking for the second 
transition towards background geometries, where only 
the shortest edge length j = 0 is allowed: Both for 
k = 4: and k = 8 this transition occurs always first for 
the (0, jmax) model followed by the (0, jmax) model, 
which coincides for fc = 4 with the J = 2 model. For 
k = 8 even all J > 0 models show transitions at larger 
Po than the other two models. 

Before discussing the area coupling, let us briefly sum¬ 
marize the results for the length coupling thus far: For 
all topological backgrounds, we observe two phase tran¬ 
sitions, one purely for the Ising model, the other in the 
geometry towards configurations, which only allow for 
j = 0, induced by the suppression of spins j > 0 by 
this particular coupling of the Ising model to the back¬ 
ground. The actual position of these phase transitions 
is primarily determined by the number of configurations 
with non-vanishing spins, not the background that allow 
for the largest spins (due to the coupling rules the quan¬ 
tum group). However, unless /3o is large enough such 
that spins j > 0 are strongly suppressed, we observe that 
the system quickly flows back to the topological hxed 
point and the actual flow of the tensor only involves the 


Ising degrees of freedom. In that sense, the Ising model 
is sensitive to the background, yet both systems are not 
strongly coupled. 


B. ‘Area coupling’ — Results 

In contrast to the ‘length’ coupling, the ‘area’ coupling 
is more non-local in nature, as we have explained in sec¬ 
tion IIC 2, as it rather prescribes a coupling with respect 
to the geometry. Still it includes the idea that the inter¬ 
action between two spins are weakened if their distance 
is increased, but it relates this distance to the areas of 
the squares this edge is shared by. Roughly speaking 

P ~ Po —-J where we simply give the area of a 

square by multiplying its edge lengths. 

As we have already stressed in section IIC 2, the mech¬ 
anism to relate the edge length to the area in the defi¬ 
nition of the modified coupling constant could prevent 
the suppression of spins j > 0 for large /3o, however, due 
to the non-local nature of this coupling, a statement for 
an edge alone is not conclusive. Indeed, this coupling 
shifts the perspective rather towards the ‘regularity’ of 
the configuration as it assigns generically larger weights 
to vertices for which the four spins je do not agree - the 
largest weight is assigned to spins alternating between 
j = 0 and j = jmax, the smallest to equilateral configu¬ 
rations for any j. 

All spins J < jmax allowed: We start the discussion 
of the results again for the background models labelled 
by a spin J, i.e. those models that allow all spins j < J. 
In table V we have summarized the results: 

The results are qualitatively very different to the 
length coupling case (compare to tables I and II). Most 
notably, we only find the phase transition for the Ising 
model from the disordered phase to the ordered one, the 
background model does not change. 

There exist geometric transitions for very larger /Iq, 
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Transition in d“interval 

fc = 4 

fc = 5 

k = 6 

k = 7 

fc = 8 

J = 0 

[0.388,0.389] 

[0.388,0.389] 

[0.388,0.389] 

[0.388,0.389] 

[0.388,0.389] 

J= 1 

[0.377,0.378] 

[0.377,0.378] 

[0.378,0.379] 

[0.378,0.379] 

[0.378,0.379] 

J = 2 

[0.373,0.374] 

[0.374,0.375] 

[0.375,0.376] 

[0.375,0.376] 

[0.376, 0.377[ 

J = 3 

n.a. 

n.a. 

[0.373, 0.374] 

[0.374,0.375] 

[0.374, 0.375[ 

J = 4 

n.a. 

n.a. 

n.a. 

n.a. 

[0.373, 0.374] 


TABLE V. Interval for parameter j3 marking the phase transition from the ‘disordered’ to the ‘ordered’ phase (on the geometry 
with all spins je < J allowed) in the ‘area’ coupling. 


e.g. po > 3 for the J = jmax models or even /3o > 10 
for the J = 1 models. This is due to the divergence 
of the edge weights as /3o —t oo (if compared to all 
j = 0), greatly shifting the model away from the topo¬ 
logical background fixed point. One can identify some 
patterns, however the interpretation is less clear than in 
the length coupling and rather indicates that a modifica¬ 
tion of the coupling for large /3o is necessary or /3o should 
be restricted to smaller values. 

The second significant change to the length coupling 
is the decrease of the phase transition temperature as 
the parameter J is increased. Recall that increasing J 
allows for more background configurations, in particular 
larger spins / edge lengths get excited. Among the added 
configurations are a few equilateral ones, but these are 
clearly outnumbered by irregular configurations, i.e. con¬ 
figurations in which not all edge lengths are equal. E.g. 
consider the case fc = 4 from J = 1 to J = 2: the lat¬ 
ter adds three allowed couplings of two spins to the spin 
j = 2, namely (2,0), (1,1) and (0,2). On the level of 
one 4-valent tensor, this adds 3 non-equilateral colouring 
((2,2,0,0), (2,2,1,1) and (2, 0,1,1)) plus perturbations 
in contrast to just one new equilateral configurations^^. 
Since these irregular configurations come with a larger 
weight than the equilateral one, the Ising spins effectively 
see a stronger coupling such that the transition occurs at 
smaller /3o. Yet this shift in the phase transition parame¬ 
ter is much smaller than in the length coupling, the most 
significant jump is the one from J = 0 to J = 1. 

Also, if we keep J fixed and instead increase k, the 
phase transition parameter changes only slightly, this 
is due to two effects. One effect applies to the models 
J = if / = Jmax = I as /c is increased by 1. As fc is 
increased j' is no longer of quantum dimension one and (a 
few) new couplings are allowed, e.g. from fc = 4 —)■ /c = 5 
the couplings (1,2), (2,1) —>• 2 become possible. Thus 
again more background configurations are allowed, how¬ 
ever the added ones are less ‘irregular’ and thus average 
out the most irregular ones further, such that the Ising 
coupling is effectively slightly weaker. In the other cases, 
this mechanism is not present, however note that the am¬ 
plitudes of the background model also change as k gets 
increased. 


Note that the Ising weights are not equal for all perturbations. 


Only j = 0 and j = allowed: 

We continue our presentation of the results for the 
background model, in which only j = 0 and jmax are 
excited for even k, i.e. only representations with quan¬ 
tum dimension one. This is particularly interesting for 
the area coupling, since this background only supports 
the smallest and the largest coupling, such that we can 
expect to see an emphasis on very irregular background 
configurations. The results for the phase transition of the 
Ising model are summarized in table VI: 

As expected, only allowing configurations with j = 0 
and j = jmax affects the phase transition parameter 
greatly: The parameter always occurs for a smaller pa¬ 
rameter than for any of the j < J models and decreases 
significantly as k, and thus jmax are increased. This is 
straightforward to understand, since the coupling rules 
and thus the allowed geometric configurations do not 
change, but due to the growth of jmax these configura¬ 
tions are more irregular such that the Ising spins experi¬ 
ence a stronger coupling. Also note that for this model no 
geometric transitions occur, since the excited representa¬ 
tions are of quantum dimension dj = I and cannot cou¬ 
ple among themselves to any other representation. More¬ 
over, in contrast to the length coupling, the area coupling 
does not suppress large spins j as /3o is increased. 

Only j — 0, j = and j = jmax allowed: Even¬ 
tually let us study the last background model, defined 
only if fc is a multiple of four, for which only j = 0, 
j = = I and j = jmax are excited. We can consider 

this model to lie ‘in between’ the extreme case of only al¬ 
lowing j = 0 and jmax and the j < J models. Again note 
that the model (0, jmax) agrees up to signs with the 
model J = 2 for k = A. The results are summarized in 
table VII: 

The results agree with our expectations: for A: = 4 the 
results are consistent / identical to the J = 2 model, for 
k = 8 we find the transition at lower /3q than for any 
j < J model, but clearly higher than for the (0,jmax) 
model. We can interpret this as follows: The addition 
of one intermediate spin, here with respect to the 
(Oj jmax) model adds many configurations that are in be¬ 
tween the equilateral and the most irregular configura¬ 
tions (alternating between j = 0 and j = jmax)- Yet 
there are only very few of these extreme configurations, 
such that they get averaged out by the majority of in¬ 
termediate configurations. This effect is enhanced once 
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Transition in dHnterval 

k = A 

k = 6 

k = 8 

Ising ‘disordered’ —^ ‘ordered’ 

[0.36,0.361] 

[0.348, 0.349] 

[0.338, 0.339] 


TABLE VI. Results for the background with only j — 0 and j = jmax allowed for even k together with the Ising model in the 
‘area’ coupling. We hnd only one transition, purely for the Ising model from ‘disordered’ to ‘ordered’ as Pq is increased. 


Transition in /3-intervaI 

fc = 4 

fc = 8 

Ising ‘disordered’ —> ‘ordered’ 

[0.373,0.374] 

[0.367,0.368] 


TABLE VII. Results for the background with only j ~ 0, j = and j = jmax allowed for fc = 4Z, / € N, together with the 
Ising model in the ‘area’ coupling. We find only one transition, just for the Ising model from ‘disordered’ to ‘ordered’ as Po is 
increased. 


more representations are allowed, as in the j < J mod¬ 
els. Again, no geometric transitions occur: even though 
the excitation of other spins is allowed, the model always 
flows back to the (j = 0, j = j = background. 

Before we continue with the discussion of the Ising 
model coupled to a superposition of the background in 
section V, let us summarize the results for the area 
coupling: As anticipated from our discussion in section 
IIC 2, this coupling shifts the focus away from just the 
length of edges to the (ir)regularity of the lattice, i.e. the 
relative size of neighbouring edges meeting at the same 
vertex, where the Ising spins effectively see a stronger 
coupling for irregular configurations than for equilateral 
ones. Considering the dynamics, the most irregular con¬ 
figurations push the phase transition of the Ising model 
significantly towards smaller Pq, as we see from a direct 
comparison of the (0, jmax) model to the standard Ising 
model (J = 0). However, both the equilateral and the 
most irregular configurations are only a small subset of 
configurations as soon as more representations je are al¬ 
lowed. Therefore, the averaging over less irregular back¬ 
grounds, like in the (j = 0, j = j = jmax) model 
and more profoundly in all j < J models (for J > 0), 
shifts the phase transition closer to the equilateral case 
and appears to be rather stable as the quantum group 
level k is increased. Again, as for the length coupling, the 
background is very stable, i.e. the system quickly flows 
back to the topological fixed point; the system mainly 
flows in the Ising degrees of freedom. Thus, while the 
Ising model is sensitive to the different background con¬ 
figurations, both systems are not strongly coupled. This 
will change in section V, where we also superimpose back¬ 
ground fixed point intertwiners, allowing these degrees of 
freedom to flow under renormalization as well. 


All fixed point intertwiners have the same dependence on mag¬ 
netic indices. Thus the transformation to the recoupling basis is 


V. THE ISING MODEL COUPLED TO A 
DYNAMICAL BACKGROUND 


After the examination of the Ising model coupled to 
a topological background, we generalize the background 
model by superposing (two) intertwiner models, which 
just amounts of superposing two different choices of pa¬ 
rameters a{{je}) discussed in section II A. As a result, the 
background model alone is no longer a topological theory 
and will flow under renormalization. In this section, we 
therefore address the following points: First, we are inter¬ 
ested to see, how the presence of the Ising model affects 
the flow of the background model, e.g. on the position 
of the phase transition between the two topological fixed 
points as Pq is changed or on the appearance of back¬ 
ground model different from both initial ones. Second, 
in section IV we have observed significant effects on the 
phase transition of the Ising model due to the presence 
of the background, such that one may wonder how it re¬ 
acts to a superposition of backgrounds. As a last point, 
the region in parameter space in which both the Ising 
model and the background model are close to a transi¬ 
tion is interesting to investigate, yet it is very likely that 
an algorithm keeping more singular values than the one 
prescribed in section III AI has to be used. 

Let us briefly outline the changes to the initial tensor 
due to the superposition of intertwiner models: The ini¬ 
tial tensor T written in the recoupling basis, see eq. (4.2), 
is simply modified by summing over the factors a({je}) 
associated to the topological background fixed point^^, 
to differentiate them we added a superscript: 


unaffected by the superposition. 
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= ^ j4,j5) + (1 - «) j2, J5)a^^^(j3, j4, js)) 

\/wfei(ji; {je})wfc2(j2; {jeDw/caOs; {je})wfejj 4 ; {je}) , 


X 


(5.1) 


where a S [0,1] is the superposition parameter. In prin¬ 
ciple one can also discuss superpositions of more models, 
we restrict ourselves to two models to present the results 
in a phase diagram for a and /3o. Note that in order 
to obtain the initial 3-valent tensor, one has to perform 
a singular value decomposition on (5.1); the sum over 
factors a{{je}) does not factorize over 3-valent vertices. 

Recalling the background models introduced in sec¬ 
tion II A, there are many possibilities to choose from; too 
many to cover in this article. Hence we will restrict the 
discussion to the superposition of two models from the 
j < J class, namely J = 1 and J = 3 for fc = 4. This is 
an interesting choice, in particular since the model J = 3 
allows the spins j = 2 and j = 3 in contrast to the 
J = 1 model. One can expect a possible transition to 
the model J = 2 between them, even without the Ising 
model. Moreover, the models J > 0 are preferable to 
e.g. the (0, jmax) or the J = 0 model, since the latter 
have very restrictive coupling rules, strongly suppressing 
the excitation of other spins. Of course, one can also en¬ 
hance the study by going to larger levels k of the quantum 
group, however this significantly increases the numerical 
effort. 

Before discussing the results for the two different cou¬ 
plings, let us briefly focus on the impact of the superpo¬ 
sition of the J = 1 and the J = 3 model at /3o = 0, i.e. 
without the Ising model, to get an idea of the geometries 
preferred by these geometries. To do so, let us study 
the singular values associated to the spins j = 1, 2, 3 
for varying a obtained from the first singular value de¬ 
composition splitting the 4-valent tensors (see eq. (5.1)) 
into 3-valent ones for the triangular algorithm in fig. 5. 
As expected, we observe the quick decline of the singu¬ 
lar values associated to the spins j = 2 and j = 3 as 
the system moves away from the J = 3 to the J = 1 
model. The fate of the spin j = 1 is more interesting: 
even though it is excited in both models, it goes through 
a minimum close to a = 0.6, which actually is very close 
to the phase transition between these two models. In 
fact this transition is not direct but rather interrupted 
by a tiny phase of the J = 2 model spanning the interval 
a S [0.60567,0.60568]. The interpretation of this data 
is rather straightforward: As a grows, geometries con¬ 
taining spins j = 2 and j = 3 become less probable and 
eventually vanish as a = 1. Configurations containing 
j = 1, which is excited in both superimposed models, 
become slightly disfavoured for intermediate values of a, 
which conversely means an emphasis on the only other 
remaining spin j = 0. Thus close to the phase transi¬ 
tion, we can expect that configurations containing j = 0 


SV(i) 

SV(j=0) 

1 . 


.. . 


^■ = 1 


*• i = 2, 3 


FIG. 5. Singular values for j = 1, 2, 3 in the first splitting 
of the 4-valent tensors, normalized with respect to the one 
for J = 0, for the superposition of the models J = 1 (a = 1) 
and J = 3 (a = 0). The values for j — 2 and j — 3, which 
have identical singular values, fall off quickly as a is increased, 
j = 1 has a minimum around a = 0.6, which is actually close 
to the phase transition (s) from J = 3toJ = 2toJ = l (for 
growing a). 


are favoured over configurations without j = 0, however 
this effect is too weak to cause the model to flow to the 
phase, in which only j = 0 is allowed. We will see that 
the Ising model, more precisely its phase transition pa¬ 
rameter, is sensitive to these changes in the background 
in both coupling schemes. 

In sections V A and V B we present and qualitatively 
discuss the phase diagrams for the length and the area 
coupling respectively. 


A. Length coupling 

The phase diagram for the length coupling can be 
found in fig. 6. Let us first start with some basic ob¬ 
servations: The vertical lines a = 0 and a = 1 are the 
pure intertwiner models for J = 3 and J = 1 respectively, 
the results are consistent with section IV. The horizontal 
line /3o = 0 represents the pure background model, since 
the Ising model assigns constant amplitudes to all back¬ 
ground configurations. There we also observe the two 
geometric phase transitions mentioned above (which are 
difficult to resolve). 

In general we observe four ‘main’ phase transitions: 
the two horizontal lines indicate the transitions induced 
via the Ising model (coupling), the lower one for the Ising 
model transition from disordered to ordered, the upper 
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FIG. 6. Phase diagram for fc = 6, superimposing the mod¬ 
els J = 1 (a = 1) and J = 3 (q = 0) in the length cou¬ 
pling. The black dots show the actual position of the phase 
transition; the connecting lines only indicate to which phase 
transition they belong. The white areas between the transi¬ 
tions belong to the same phase and are named accordingly. 
The J — 2 phases are located between the green / blue lines, 
which is indicated by the arrows. The colours are chosen to 
according to the transition: red indicates the transition of 
the Ising model, blue the geometric transitions in the ordered 
Ising phase, green in the disordered one. Orange indicates the 
transition to the ‘just j = 0’ phase. 


one marks the transition to geometries, for which only 
j = 0 is allowed, induced by the suppression of j > 0 in 
the matter coupling. The two vertical lines, which are 
mostly indistinguishable, mark the two geometric transi¬ 
tions via the J = 2 model. Let us discuss these transi¬ 
tions separately: 

Considering the Ising transition from the disordered to 
the ordered phase, we generically observe a drop in the 
phase transition parameter as we approach the geometric 
transition, where the J = 3 line appears to merge with 
the geometric transitions before meeting the line from 
J = I. This drop in phase transition temperature implies 
that the Ising spins effectively see a stronger coupling as 
one moves away from the topological fixed points, where 
this is far more striking compared to the J = 3 perspec¬ 
tive. Considering the properties of the length coupling, 
the increase in the effective coupling for the Ising spins 
implies an emphasis on smaller spins in the background 
configurations. As we have concluded above from fig. 
5, this is exactly what happens: Seen from a = 0, i.e. 
the J = 3 model, the spins j = 2 and j = 3, which 
weaken the interaction between neighbouring Ising spins, 
are significantly suppressed as a is increased explaining 
the substantial drop in the phase transition parameter. 
Moreover, even the spin j = 1, which is excited in both 
superimposed models, gets disfavoured close to the phase 
transition, such that also the effective coupling is weaker 
than in the J = 1 model. 

The story is similar for the second transition induced 
by the Ising model, namely the one to ‘just j = 0’ ge¬ 
ometries. Again we observe a drop in the phase transi¬ 


tion parameter from both sides close to the (significantly 
shifted) geometric transition. Recalling that this transi¬ 
tion occurs due to the suppression of all spins j > 0 by 
the Ising edge weights, it is clear that this suppression 
occurs for smaller /3q if configurations allowing just j = 0 
become emphasized by the geometric superposition. 

Eventually, let us discuss the geometric transition: as 
mentioned above, there are actually two geometric tran¬ 
sitions, enclosing the emergent J = 2 configuration. This 
pair of transitions starts out vertically from the x-axis, 
however quickly bends toward smaller a, i.e. J = 3, 
where it also becomes broader after crossing the phase 
transition line for the Ising model. Again this is straight¬ 
forward to explain from the length coupling of the Ising 
model: As P{j) falls of as y, larger spins j get quickly 
suppressed. The spin that is suppressed the most is 
j = 3, then j = 2, etc. Thus close to the geometric 
phase transition this suppression favours a flow towards 
the J = 2 background over the J = 3 background and 
analogues for the J = I and J = 2 background. If the 
system is then moved to smaller a, j = 3 in turn be¬ 
comes emphasized by the background such that the sup¬ 
pression from the Ising weights must be stronger for the 
J = 2/J = I phase to appear, which explains the ad¬ 
vancing of both the J = 2 and J = I phase to smaller a. 
At some point, the suppression is too strong, such that 
all spins j > 0 are suppressed and the geometric transi¬ 
tion ‘merges’ into the transition towards geometries, for 
which only j = 0 is allowed. 

To conclude both the Ising and the background model 
react very sensitively to one another. From the effective 
coupling (and thus the position of the phase transition) 
seen by the Ising model, one can deduce which geometric 
configurations are preferred at a specific superposition 
of background intertwiner models. Here the superposi¬ 
tion results in an emphasis towards edges carrying spin 
j = 0, such that the effective coupling for neighbour¬ 
ing Ising spins is stronger. Conversely, the geometry is 
greatly affected by the matter coupling as well: if the 
Ising coupling is strong enough, a geometric transition 
to a different background model is induced. Admittedly 
this requires the background to be in a superposition, 
otherwise the deviation away from this fixed point caused 
by the Ising model is not sufficient. Also note that we 
did not find any phases beyond the ‘product phases’ of 
the combined model. Whether more phases exist, e.g. 
close to both phase transitions, can probably only be an¬ 
swered by using a tensor network algorithm keeping more 
singular values in each iteration, yet close to the phase 
transition an infinite bond dimension is required. 

In section VB we discuss the significantly different 
phase diagram for the area coupling. 


B. Area coupling 

The phase diagram for the area coupling can be found 
in fig. 7. As expected, we observe a very different phase 
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J = 1, ordered 

J = 3,ordered 

4 —J = 2, ordered 

J = 3, disordered 


J = 2, disordered 

J = 1, disordered 



0.606 0.608 0.610 0.612 0.614 

FIG. 7. Phase diagram for fc = 6, superimposing the models 
J = 1 (a = 1) and J = 3 (a = 0) in the area coupling. 
The plotstyle, in particular the colouring of the transitions, is 
similar to the one for the length coupling. 


diagram in comparison to the length coupling. First of 
all, the Ising model transitions for a = 0 (J = 3) and 
a = 1 (J = 1) are consistent with the results in section 
IV. Again we observe two geometric phase transitions 
for /3o = 0, in the same small interval as for the length 
coupling. 

The horizontal line marking the phase transition of 
the Ising model from the disordered to the ordered phase 
(as /3o is increased) is indeed almost horizontal, it only 
slightly decreases (from both sides) to roughly /3o « 0.371 
as the geometric phase transition is approached. Thus 
the Ising spins again see an effective coupling that is 
slighly stronger as the geometric transition is approached; 
as the phase transition occurs for larger /3o for a = 1 
(compared to a = 0), we can also conclude that this 
effect is stronger in the regime where the J = 1 back¬ 
ground dominates. Let us analyse this behaviour in 
steps: Recalling the discussion for the superposition of 
backgrounds we already know that the superposition of 
intertwiners affects the amplitudes such that the trivial 
representation j = 0 is emphasized, while all others get 
suppressed, albeit only slightly in the case of j = 1. As 
we have already discussed before the area coupling is less 
sensitive to the excited edge lengths alone, but rather 
how irregular the geometry is as seen from a single ver¬ 


tex. Indeed, the strongest coupling of the Ising spins 
(for a single vertex) occurs if the edges meeting at this 
vertex alternate between minimal (j = 0) and maximal 
(j = jmax) length. Aware of this mechanism, the em¬ 
phasis of J = 0 close to the geometric transition explains 
the slightly stronger effective Ising coupling as follows: 
If j = 0 is slightly favoured, while all also spins j > 0 
remain excited, this puts an emphasis on more irregu¬ 
lar configurations. Additionally, from the perspective of 
the J = I model, the superposition excites spins _) > I as 
well, such that more irregular configurations are possible, 
resulting in the slightly larger decline as the geometric 
transition is approached from a > 1. 

Similar to the transition of the Ising model, the two 
geometric transitions, i.e. from from J = 3 to J = 2 to 
J = I as a is increased, is also less affected by the cou¬ 
pling to the Ising model than in the length coupling. The 
most crucial difference is the ‘bending’ of the transition 
towards larger a instead of smaller a. Again this is due 
to the different nature of the area coupling as it is rather 
sensitive to the (ir)regularity of geometries than just the 
excited edge lengths. As /3o is increased, the irregular 
geometries are given more weight than the regular ones 
from which one can conclude a favouring of larger spins, 
since these allow for more irregular geometries, in par¬ 
ticular with respect to j = 0. Thus the transition to 
models with J > 1 occurs at larger a, yet this effect is 
much weaker than in the length coupling. Also, we do 
not observe a broadening of the J = 2 phase similar to 
the situation in the length coupling. 

Again, we would like to mention that additional ge¬ 
ometric transitions occur for very large j3o due to the 
divergence of the Ising weights, which occur for J = 3 
at roughly /3o > 3 and for J = 1 at roughly /3o > 10. 
We will not cover these in detail in this work because 
of their complicated structure and since they occur in 
a regime far away from the interesting dynamics of the 
Ising model. To us this rather indicates that further mod¬ 
ifications to the matter couplings must be done if these 
are to be avoided or /3o should be restricted to smaller 
values. 

Let us some up the results for the area coupling. Again 
we observe that both the Ising model and the geomet¬ 
ric background are sensitive to each other, however to a 
much smaller degree than for the length coupling. The 
effect on the position of the Ising phase transition is com¬ 
paratively small and can be well explained by the changes 
in the background geometry induced by the superposition 
of fixed point intertwiners, namely the slight preference 
of irregular geometries due to an emphasis of configura¬ 
tions containing edges with spin j = 0. In contrast to the 
length coupling, the geometric transitions shift towards 
larger a, i.e. towards the J = 1 model, as /Iq grows. 
This is caused since for larger /3o more weight is assigned 
to irregular configurations, which prefers the appearance 
of spins j > I, even though these are disfavoured by the 
background model. Again the coupling to the Ising model 
induces geometric transitions, here from J = I to J = 2 
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and then J = 3. As already mentioned for the length 
coupling, a more thorough study might be worthwhile 
at the intersection of both transitions, yet this requires 
an algorithm capable of efficiently studying systems at 
criticality, which means obtaining good results, e.g. for 
critical exponents, at a finite bond dimension. 


VI. SUMMARY AND DISCUSSION 

In this article we have studied a 2D toy model for 
coupling pure lattice Yang-Mills theory to a spin foam 
model, where the gravitational and the gauge theory are 
replaced by simpler models, which do preserve parts of 
the dynamical structure of the full theories. For the 
gravitational part, we have chosen topological intertwiner 
models defined in [59], which allow for a notion of inter¬ 
twiner degrees of freedom similar to the spin foam case, 
where these occur due to the imposition of simplicity 
constraints. These intertwiner models are defined for 
the quantum group SU(2)fc, which comes with a natu¬ 
ral cut-off on the representation labels [60, 61]. More¬ 
over it is conjectured that spin foam models defined for 
quantum groups incorporate a non-vanishing cosmologi¬ 
cal constant [78-90]. For the matter system, we have con¬ 
sidered spin systems, more precisely the Z 2 Ising model, 
for which the local gauge symmetry of lattice Yang-Mills 
is replaced by a global symmetry. Again it is expected 
[28] that these models share statistical properties with 
their 4D lattice gauge theory relatives. 

Due to the structure of the two systems it is straight¬ 
forward to coarse grain them via tensor network renor¬ 
malization [57, 58]. To do so the partition function of 
the system is rewritten as a contraction of a tensor net¬ 
work. By performing local manipulations of the tensors, 
in particular defining and identifying the most relevant 
coarse degrees of freedom via a singular value decomposi¬ 
tion, the partition function is evaluated / approximated 
in steps until eventually the procedure reaches a fixed 
point. In context of this coarse graining procedure, we 
have briefly revisited the discussion in [24] whether the 
matter degrees of freedom should be placed on the same 
lattice as the spin foam or rather the dual lattice as it is 
done e.g. in [22]. For the toy model under discussion, this 
choice has rather technical than physical consequences: 
First, the Ising model is self dual [67], i.e. both choices 
are related by a duality transformation, which maps the 
strong coupling sector to the weak coupling sector and 
vice versa. However, in order to coarse grain these mod¬ 
els, they have to be cast into tensors, in essence local am¬ 
plitudes that are associated to the vertices of the lattice. 
As it turns out, if the Ising model is placed on the dual 
lattice, it is more suitable to keep the original Ising spin 
variables; in the process they have to be coarse grained 
like decorated tensor networks [50] , where the Ising spins 
‘decorate’ the gravitational part. Conversely, if the Ising 
model ‘lives’ on the same lattice, it is preferable to group 
Fourier transform the Ising variables. These insights for 


coarse graining coupled systems might be helpful, e.g. 
to examine the influence of the choice of discretisation, 
once one tackles more complicated theories, in particu¬ 
lar if these theories are not self dual or the dual theories 
are not known. Nevertheless, we agree with [24] that it 
is preferable to define both matter and gravitational de¬ 
grees of freedom on the same discretisation, such that one 
can establish contact to the states and Hilbert spaces of 
the associated canonical theory. 

Undoubtedly, the crucial part of this article is the mat¬ 
ter coupling between the background system and the 
Ising model: The background intertwiner models endow a 
geometric interpretation onto the lattice, its SU( 2 )fe rep¬ 
resentation labels je on the edges assign a length je + 5 to 
these edges. Intuitively, the interaction among Ising spins 
should be sensitive to the distance between them, spins 
further apart should have a weaker interaction than close 
spins. Therefore we have introduced a coupling between 
Ising spins and the background by choosing the Ising cou¬ 
pling constant /? to be a function of the background spins 
j, actually similar to [22, 23]. In fact, it is also related 
to the idea expressed in [47] to assign local coupling con¬ 
stants (to parts of the discretisation) for lattice gauge 
theories defined on irregular lattices. Interestingly, the 
related construction of the Ising model on 2D causal dy¬ 
namical triangulations [55, 56] is somewhat orthogonal to 
our idea: Instead of fixing the discretisation, in particu¬ 
lar the combinatorics, and varying the edge lengths, one 
rather fixes the edge lengths and sums over all possible 
triangulations. Thus, the distance between neighbouring 
Ising spins is always the same, thus no modification of 
the coupling constant /?, but the number of neighbouring 
Ising spins can vary. 

In this article, we have constructed two modifications 
of the coupling constant, one for which /3(j) ^ j and 
one for which the distance is related to the (square 
root of the) area that is partially spanned by this edge 

/3(j) ^ While the first is a local interac¬ 

tion, assigning the highest weight to the shortest edge 
length, the latter is more non-local by relating the edge 
length to the overall geometry and is analogous to the 
modification introduced in [22] for 4D Yang-Mills cou¬ 
pled to the Barrett-Crane spin foam model [ 68 ]. There 
it was chosen such that the Wilson action scales as the 
volume of the 4-simplex it is defined on. Admittedly, 
these modifications are heuristic and appear ad hoc; one 
can certainly construct more realistic couplings, e.g. by 
considering Laplacians on discretisations, and it would 
be preferable to derive such a coupling from the classical 
theory in 4D. However, neither is this the main purpose 
of this article, nor does this avoid choosing a discretisa¬ 
tion, which eventually results in lattice imperfections. It 
is rather to explore whether a simple modification in the 
matter part, which couples the two systems, can result in 
an interesting and reasonable dynamics for the coupled 
system or whether a different mechanism is necessary. 

Indeed, the chosen couplings have a significant effect on 
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the dynamics as they also emphasize different geometries 
of the background system. The straightforward length 
coupling simply considers edge lengths, the longer an 
edge the weaker the interaction between Ising spins is. As 
a result, the Ising model sees effectively a (much) weaker 
coupling the more edge lengths j > 0 (and thus also con¬ 
figurations of the background with j > 0) are allowed, 
such that the phase transition is clearly shifted towards 
larger /3o- However, this coupling has a direct effect on 
the background geometry as for /3o 1 spins j > 0 get 

strongly suppressed causing a phase transition to geome¬ 
tries with just spins j = 0 allowed. Moreover, the Ising 
model is very sensitive, in terms of the position of its 
phase transition, to changes in the background model. 
In particular the emphasis on the smallest possible edge 
length suggests two possibilities: It can imply that the 
coupling /?o should be restricted to smaller values or one 
should introduce a modification to the coupling. 

In contrast to that the area coupling is not very sensi¬ 
tive to edge lengths alone; it rather favours specific con¬ 
figurations due to its non-local nature. Interestingly, it 
assigns the weakest weights to regular, equilateral config¬ 
urations (independent on the edge lengths as it is invari¬ 
ant under a global rescaling of all edge lengths) and the 
strongest to very irregular configurations, e.g. if around 
one vertex the edge lengths alternate between the min¬ 
imal and maximal edge lengths. As a consequence, the 
Ising spins effectively see stronger couplings if spins j > 0 
are allowed, thus shifting the phase transition to smaller 
j8o, yet this shift is considerably smaller in comparison 
to the length coupling and far less sensitive to changes 
of the background model. E.g. when considering a su¬ 
perposition of background models in section V, the line 
marking the Ising phase transition is almost horizontal 
with only slight deviations consistent with the changes 
in the background geometry. In fact this is encouraging, 
since we do observe that this coupling allows for matter 
sensitive to the geometry it is defined on, yet without re¬ 
sulting in significant deviations away from the equilateral 
case. 

Fortunately, also the gravitational background model 
is sensitive to the matter defined on it, albeit only signif¬ 
icantly if it is not close to a topological fixed point. If the 
latter is the case, the deviations away from it caused by 
the Ising model are not strong enough, the background 
quickly flows back to the initial fixed point intertwiner 
and decouples from the Ising model, which remains sen¬ 
sitive to the background geometry. However, if we con¬ 
sider the superposition of two background fixed points, 
the interaction with the Ising model can trigger phase 
transitions in the geometry. In fact, depending on the 
geometries preferred by the coupling scheme, for growing 
/3q the geometric transition is shifted away from its posi¬ 
tion in parameter space without the Ising model. Again 
this effect is far more profound in the length coupling 
than in the area coupling. Nevertheless, we have shown 
that by introducing a simple modification of the Ising 
model coupling constant, which couples it to a topolog¬ 


ical background, one can implement a dynamics, which 
exhibits regions (in parameter space) of strong coupling 
between the two systems, namely close to the phase tran¬ 
sitions of both systems, and regions in which they are 
only slightly coupled, e.g. on the fixed points of the back¬ 
ground model. Hence, it might be worthwhile to study 
a similar coupling mechanism for spin foams and lattice 
gauge theories in order to potentially identify similar re¬ 
gions as well. A region in which both systems are only 
weakly coupled might allow us to establish contact to 
lattice gauge theories on a fixed (flat) background, which 
could indirectly teach us more about the geometry arising 
from spin foam models. 

As these coupling mechanisms of the Ising model to 
the background are the central piece of this article, let us 
also discuss some of the choices made and their possible 
consequences. In both cases, we have regularized the 
j = 0 case by defining the j + | as the length of an edge. 
We have motivated this choice with the asymptotics of 
(3D) spin foam models [69, 70], but admittedly this is ad 
hoc and one can straightforwardly generalize this to j -I- e 
instead, where e ^ 0. 

For the length coupling, e > 0, i.e. essentially the case 
studied here, implies that /3(j) is always positive and the 
largest coupling constant is assigned to j — 0, and thus 
also the largest Ising edge weight. The size of e then de¬ 
termines how much j = 0 is favoured, but the qualitative 
behaviour is more or less the same. However, if e < 0, 
/3(j) is negative for all j < jej, such that for some edge 
lengths, the interaction for j < jej is antiferromagnetic 
and ferromagnetic for j > jej. One can expect signifi¬ 
cant changes for the dynamics of such a model, if one 
tunes e such that small edge lengths give antiferromag¬ 
netic and large edge lengths ferromagnetic interactions 
between the Ising spins. One can also choose e such that 
the interaction is antiferromagnetic for all j, where |/3(j)| 
is then largest for jmax- Thus it might be possible to 
study other interesting phenomena, possibly a transition 
from antiferromagnetic to ferromagnetic or the modified 
interaction could cure the unintended emphasis on j = 0 
configurations. However one should not literally interpret 
j -I- e as the length, but rather an effective description, in 
which the system rather avoids shortest distances, which 
one could associate with highly energetic configurations. 

Concerning the area coupling, a modification to e < 0 
can also have interesting effects, which are however not 
as straightforward to predict, as /3(j) can become com¬ 
plex (due to the squareroot). Another option can be to 
invert the modification, i.e. normalize the length by the 
square root of the area instead. However, it is not obvi¬ 
ous whether this leads to a favouring of equilateral over 
irregular configurations, as both the product of weights 
around a single vertex remain unchanged for the equilat¬ 
eral and the alternating case. This has to be examined 
more thoroughly before one can draw a conclusion. 

Before we conclude, we would like to make a few 
remarks on the coarse graining algorithm. Since we 
have used a symmetry protecting version similar to 
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those introduced in [51] for Abelian and in [53, 54] for 
non~Abelian and quantum group symmetries, we have 
been able to straightforwardly identify the phases of the 
model. The main advantage of the triangular tensor net¬ 
work algorithm, that we have used in this work, is its 
improved memory and computational cost, as we only 
save 3-valent tensors and contract one index less. This 
roughly reduces the computation time from y® to y® com¬ 
pared to the 4-valent algorithm in [58], where y is the 
bond dimension (or index range) of the tensor. However, 
the algorithm has a caveat: Even though we have worked 
at rather low accuracy, the fixed point tensors have shown 
so-called CDL structure. This CDL structure means that 
part of the tensor is of a corner double line form, its edges 
being double indices, which are pairwise identified along 
the vertex. This structure is a fixed point of many ten¬ 
sor network renormalization schemes and for fixed initial 
parameters, the final fixed point tensor is cut-off depen¬ 
dent, see also the appendix of [77] for a nice explanation. 
Thus it obstructs a proper flow to the true fixed point 
of the system. In condensed matter physics, this struc¬ 
ture is interpreted as short range entanglement, which is 
unintentionally promoted to a larger scale by the coarse 
graining procedure and can be cured by entanglement 
filtering [58, 77, 91, 92]. Due to the absence of a (back¬ 
ground) scale, the interpretation in spin foam models is 
less clear. Fortunately in the system under discussion 
here, the different phases can still be clearly identified, 
such that the results (in the approximation scheme) are 
consistent. 

Furthermore, we would like to briefly comment on the 
chosen cut-off scheme, namely one singular value per in¬ 
tertwiner channel {j,k). Even though it appears to be 
a strong simplification, our results suggest that this ap¬ 
proximation is good enough to capture many important 
aspects of the system. Therefore, if we would work with 
a higher accuracy, i.e. more singular values are taken into 
account in each iteration, we expect that the actual posi¬ 
tion of the phase transitions will change, however not the 
qualitative behaviour, e.g. which phases the model has 
or the qualitative reaction of the Ising spins interaction 
to the background. Close to both the Ising and geomet¬ 
ric transition, i.e. where both systems are more strongly 
coupled, this approximation breaks down and one has to 
work with a larger bond dimension. There we also ex¬ 
pect interesting dynamics to occur, possibly with addi¬ 
tional phases beyond the product phases of the composite 
model. Furthermore, it would be interesting to compute 
critical exponents in this regime in order to learn more 
about the phase transitions. However many tensor net¬ 
work algorithms, including the one used here, require an 
infinite bond dimension to study systems at criticality 
and are not well suited for extracting critical exponents 
[93]. The algorithm recently developed by [77] is more 
promising to fulfil this task as it successfully deals with 
CDL structure, may allow for good approximations close 
to criticality for finite bond dimension and is claimed to 
give good results for critical exponents. Therefore, it is 


more promising to examine the toy model at criticality 
with a properly adapted algorithm. 

Additionally, one may want to supplement the numer¬ 
ical analysis of the toy model by more analytical means, 
e.g. examining whether one can define the toy model to 
be discretisation independent and thus possibly fix the 
coupling of Ising spins and the background. Therefore 
one could attempt to find a model independent under 
Pachner moves, or at least perturb around the topological 
fixed points of the background and determine the Ising 
weights, such that the system is approximately discreti¬ 
sation independent. In fact we have examined these pos¬ 
sibilities, however did not find a straightforward solution 
to either of them. One rather contrived option would be 
to construct the whole amplitude, that is Ising plus back¬ 
ground, to be a topological fixed point from [59], which 
however requires the background to completely absorb 
the dependence on the Ising degrees of freedom. For the 
other option, the perturbation around a topological fixed 
point, the equations essentially force the Ising model to 
assign constant weights to all background configurations; 
the Ising model is forced to either be in the ordered or 
disordered phase. This is also reflected in the numeri¬ 
cal simulations for the topological background in section 
IV, where the system quickly flows back to these and 
only the Ising spins flow to either the ordered or disor¬ 
dered phase. In fact, this is already non-trivial to define 
the Ising model itself in a discretisation independent way 
(away from either the ordered or disordered phase), since 
one can expect it to be non-local as one can already de¬ 
duce from a simple decimation scheme and as it is also 
well-known for other interacting systems, e.g. discrete 
gravity under Pachner moves [94-97]. A possible way to 
overcome this is to generalize the concept of discretisa¬ 
tion independence by a means to relate theories defined 
on different discretisations. One way is to allow refining 
and coarse-graining of the boundary (data), in such a 
way that one can define embedding maps that relate be¬ 
tween configurations / states defined on different bound¬ 
aries. If these embedding maps are designed such that 
one can unambiguously relate different boundaries, one 
is able to identify states across them, i.e. recognize the 
same (physical) situation represented on different bound¬ 
aries. This condition is known as cylindrical consistency, 
e.g. realized for the kinematical Hilbert space of loop 
quantum gravity [9, 98, 99]. To realize this for interact¬ 
ing theories or the physical Hilbert space of loop quantum 
gravity, i.e. states annihilated by all constraints of the 
theory, is non-trivial, however see [39, 40, 100] on how 
this problem can be tackled. 

At last, we would like to raise a conceptual issue 
concerning the coupling of matter to spin foams. The 
Ising model, and also lattice Yang-Mills in the Wilson 
or heat kernel action [28], are essentially Wick-rotated 
theories that appear as exp(— 5'matter) in the partition 
function / path integral. However, spin foams, and also 
the discussed intertwiner models, are complex assigning 
~ exp(iS') to the path integral (for one orientation), also 
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for Riemannian signature. Moreover, it is well-known 
from spin foam asymptotics [101-106] that both orien¬ 
tations of a 4-simplex have to be considered, such that 
one rather finds ^ cos(S'Fiegge)) where S'uegge is the Regge 
action [107] associated to the 4-simplex. To the author’s 
best knowledge, it is not known whether it is consistent 
to couple a Wick-rotated matter theory to a non-rotated 
gravitational theory and, moreover, in which way mat¬ 
ter should couple to different orientations. There ex¬ 
ist attempts to remove the other orientation from the 
asymptotics [108-112], whereas it has been argued in 
[113] that fermionic degrees of freedom might be sensitive 
to changes in the orientation. 

In this article we have introduced and tested a toy 
model to demonstrate that coupling matter degrees of 
freedom to a dynamical (discrete) background a la spin 
foams can be straightforwardly achieved and can addi¬ 
tionally result in an interesting dynamics, where both 
systems are sensitive to one another. Moreover, we have 
identihed regions in parameter space of strong and weak 
coupling between matter and background degrees of free¬ 
dom. Therefore we suggest to generalize these ideas to 
the 4D theories, e.g. Yang-Mills theory and spin foams, 
and attempt to extract new insights for spin foam models. 
Of course these theories are very complex and a renormal¬ 
ization algorithm is (at the moment) not at hand. There¬ 
fore it will be necessary to introduce approximations and 
simplihcations, which might still allow us to study inter¬ 
esting aspects, albeit with the caveat of lacking a contin¬ 
uum limit. Nevertheless, this may allow us to tentatively 
establish connections to other approaches based on dis¬ 
cretisations, such as lattice gauge theories, or may give 
us hints on how to modify the matter coupling, the spin 
foam model or both. An application, which might be in 
reach, could be the identification of an effective (mat¬ 
ter) dynamics consistent with lattice gauge theories on 
flat spacetime, which would massively support spin foam 
models as viable quantum gravity candidates. 
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Appendix A: Quantum group basics 


In this appendix, we briefly introduce some basics on 
the quantum group SU(2)fc and the diagrammatical cal¬ 
culus, which are necessary to understand the calcula¬ 
tions, in particular in appendix B and partially in section 
III, however not essential to understand the main results 
of the paper. We use the notation and conventions in¬ 
troduced in [60] , where one can also find a more detailed 
introduction to quantum groups, see also [61]. Many of 
tools we discuss and use in this work have been originally 
developed in [54]. 

If we refer to the quantum group SU(2)/j, we are ac¬ 
tually referring to the q-deformation Uq{5u{2)) of the 
universal enveloping algebra U{su{2)) of the Lie algebra 
su(2) as in [60]. The algebra Uq{su{2)) is generated by 
three operators J ±, with commutation relations 

[Jz, J±] = ±J± 

= (Al) 

As mentioned in the main body of the paper, the finite di¬ 
mensional representation of SU(2)fc are labelled by j S | 
and can be defined on 2{j + 1 ) dimensional representa¬ 
tion spaces Vj as for SU(2). The quantum dimension dj 
of representation j is defined as the quantum number of 
the classical dimension: 


dj ■■= [23 + 1] , (A2) 


where the brackets denote quantum numbers: 


n] 



(A3) 


In this paper, the deformation parameter g is a root of 
unity, with q = exp(p^^i). A: G N is called the level 
of the quantum group SU(2)fe. Quantum numbers are 
periodic 
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with zeros at n = 0 and n = k + 2. Thus / = f with 
dk /2 = ^ is the ‘last’ representation with a strictly posi¬ 
tive quantum dimension. Representations j = 0, |,..., | 
are called admissible, representations j > | are of so- 
called quantum trace zero. 

The tensor product of two representations via 

the co-product A. The action of the SU(2)fc algebra on 
Vji 0 Vj 2 is defined as 


A(J±) = q 0 J± - 1 - J± 0 

A( J^) = 1 0 -t- Jz 0 1 . (A5) 


The tensor product Vj^ 0 Vj^ can be decomposed into 
a direct sum of irreducible representations plus a part 
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consisting of trace zero representations (which are mod- 
ded out). With an orthogonal basis \j,m) in the repre¬ 
sentation spaces, the decomposition is given by Clebsch- 
Gordan coefficients 

\j,m) = ® ^ 2 ,^ 2 ) . (A 6 ) 

mi,m 2 

If one couples three admissible representations jj, Jk 
and jL in this way, several conditions have to be satisfied 
for the Clebsch-Gordan coefficients to be non-vanishing: 

ji + jK > Jl for permutations {J,K,L} oi {1, 2, 3} , 
ji + J 2 + J 3 = 0 mod 1 , 

ji + j2 + js ^ k . (A7) 

The last condition in (A7) is special to the quantum de¬ 
formed case at root of unity and indicates that Vj, 0 1^2 
can include trace zero parts, which can be modded out 
[61]. However, some equations (for instance the defini¬ 
tion of the [ 6 j] symbol) are only valid up to trace zero 
parts [61]. 

In particular we have the completeness relation 


E 


(^jlj2j3 

Q^mim2m^ 


m3, ^3 admiss. 


rijij2j3 

y mj^m2m3 


= , ,(A8) 

mim 2 ,m-^m^^ > 


where projects out the trace zero part in 

Vj, 0 Vj 2 - The orthogonality relation for the Clebsch- 
Gordan coefficients is given as 


E nhhi fihhi' 

q^mim2m q^mim2m’ 

mi ,m2 




where dj,j 2 j = 1 if the coupling conditions (A7) are sat¬ 
isfied and vanishing otherwise. 


They are interpreted as a map Vj, 0 Vj^ —> Vj^, sym¬ 
bolizing how the spins ji and j 2 (with their respective 
magnetic indices) couple to J 3 . We have already intro¬ 
duced their graphical representation in equation ( 2 . 1 ): 

h 

■■= • (AlO) 

jl J2 

A particular version of this Clebsch-Gordan coefficient 
will be important later on: If we choose ji = J 2 = j and 
take ja = 0, we dehne the ‘cap’ as a map: Vj®Vj ^ C, 
namely 

j 

m m' 

(All) 

From this ‘cap’ we can similarly define a ‘cup’ by requir¬ 
ing that they give the identity if we concatenate them: 



C" , (A12) 

m m 

which gives: 

m m' 

^ ={-iy+^q^6^^_rn' . (A13) 



1. Diagrammatic Calculus 

The introduction of a quantum group complicates 
some definitions known in the classical case, e.g. the 
notion of a dual, which is necessary to calculate the re¬ 
coupling basis of the tensor in appendix B and section III. 
To do so, a convenient graphical representation has been 
introduced in [54] . The quantum group requires to spec¬ 
ify a special direction, which we will take as the vertical 
direction and can be interpreted as maps from a tensor 
product of representation spaces of SU(2)fe, represented 
by incoming lines from below, to a tensor product of rep¬ 
resentation spaces, drawn as outgoing lines on top. Each 
of these lines carries a representation label j and a mag¬ 
netic index m. One basic example of such a map are the 
Glebsch-Gordan coefficients, denoted by 


Using these ‘cups’ and ‘caps’, we can construct the 
Glebsch-Gordan coefficients for the quantum group with 
inverse (here: complex conjugate) deformation parame¬ 
ter q, already given in ( 2 . 2 ), by ‘bending up’ one of the 
lower legs of the Glebsch-Gordan in (AlO). 



(A14) 

This map can hence be interpreted as mapping Vj^ —>■ 
Vji 0 Vj 2 , thus it is dual to (AlO). With a ‘cap’ we can 
‘pull down’ one of the legs again and arrive back at (AlO): 


This is not the standard Clebsch-Gordan coefficient defined in 
[60], but it is modified by the quantum dimension: = 

(-'0102 33 ( 

qCymi m2 m3 V V / 
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Concatenating these two maps, we obtain a map Vj^ 
Vj^ proportional to the identity. 


J3 js 



(A16) 

Given these graphical ingredients, several important 
identities can be derived, which we append here: 

We start with the expression giving the dependence 
of the 4-valent intertwiner (2.19) on the magnetic 
indices: 


can be used to project onto the basis (A17), from which 
also follows that the components of a 4-valent tensor in 
that basis, e.g. T in appendix B / section III, can be 
computed by contracting that tensor with (A18). More¬ 
over, we can also compute the change of basis related to 
a different 4-valent recoupling scheme: 

We are looking for the relation 


In order to find the coefficient c{j^,je) we contract the 
above expression with the diagram (AI8): 







E rji 32 is i4 is 

m2 ms ms 


J3 J4 


(A17) 


Its dual is defined by placing ‘cups’ on its bottom legs 
and ‘caps’ on its top ones (such that these do not cross). 
In terms of Clebsch-Gordan coefficients this reads: 



2j5 \ ' „m 5 njx j 2 35 nji ji 35 

/ j y g'^mi m 2 ms 1714 ms ’ 

ms 

(AI 8 ) 


In fact, if we connect the diagrams (A17) and (A18), we 
obtain the following diagram 





Ji 


from which we can deduce that the map 
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(A19) 


P({m},{TO'})(ji, j2,j3,j4) := 





J 2 J 1 


J5 


35 


h ji 



32 jl 



33 3i 


(A20) 



The diagram on the left hand side is the 6j symbol of 
SU(2)fe, given as a particular contraction of four Clebsch- 
Gordan coefficients: 




Calculating both sides of (A22), we find 


(_l)4l+42+43+44 


Ji 32 35 

33 ji J6 


(A23) 


= (_l)4l+42+43+44 ^ ^ (A24) 


from which we deduce 
j2 jl 


= Vw — 

33 Ji 


35 


36 


Jl 32 35 
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J2\ /Jl 
J3^ Aj4 


(A25) 


As a last step, let us briefly explain how the diagram 
(BI4) enters into (BI2) and (BI3) and also why (BI4) 
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can be split into two 6j symbols for the 3-valent algo¬ 
rithm in section III. 

If we consider just the magnetic indices of the interme¬ 
diate 3-valent tensors Si (equations (B8) to (Bll)) and 
connect the 3-valent vertices according to (B12) or hgure 
fig. 8, we find a diagram with two bottom legs and two 
top legs connected by a square. To compute the compo¬ 
nents of (B12) in the basis (A17), it has to be contracted 
with (A18) resulting in (B13). The diagram that we have 
to compute is thus: 



Instead of doing this directly we use the following two 
identities to simplify it first: 



Thus we can manipulate (A26) in the following way: 



(A29) 


The diagrams on the right hand side are again 6j-symbols 
defined in (A23). Eventually, we find 
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Appendix B: 4—valent tensor network 
renormalization: general idea and derivation of 
symmetry protecting algorithm 


In this appendix we will explain the 4-valent tensor 
network algorithm in more detail. As mentioned already 
in section III, the general idea is to encode the entire dy¬ 
namics of a (discrete) system, in particular the partition 
function, into a tensor network, i.e. a local contraction 
of multidimensional arrays, and evaluate it in steps via 
coarse graining. The latter is done via a local manipu¬ 
lation of the network, in which one defines new coarse 
degrees of freedom from the finer ones and introduces a 
truncation, which allows for a control on the error being 
made. Therefore, it is necessary to know the relevance of 
the coarse degrees of freedom (with respect to the other 
ones). 

Such variable transformations are computed via a sin¬ 
gular value decomposition (SVD). To explain this let us 
discuss the concrete algorithm developed in [57, 58] . Con¬ 
sider a tensor network on a 2D square lattice with iden¬ 
tical 4-valent tensors Tabcd on all vertices. Divide the 
tensors into even and odd ones, and split them accord¬ 
ing to fig. 8, where the odd ones are green and the even 
ones blue. This splitting is performed via a singular value 
decomposition acting on the following matrices: 


^{ab)\{cd) '^(^ab),{cd) 


~ 2^^(ab),i^i W }(cd),i 


* {da)\{hc) 


=:M, 


( 2 ) 

(da),{be) 

(2) C2) 




(BI) 


(B2) 


where i ^{cd) i singular vectors, Af ^ the 

singular values of the matrix ^ ^ 


{ab) denotes that the indices a,b have been combined into one 






























32 



FIG. 8. Tensor network renormalization for a regular square 
lattice. The tensor T is split into 3-valent tensors in two 
different ways by grouping two exterior legs together and 
performing an SVD. A truncation is performed on the new 
indices i by only taking the largest x singular values into ac¬ 
count. The new truncated tensors Si are combined into a new 
4-valent tensors by summing over the old indices. 


unitary matrices, i.e. the singular vectors are orthonor¬ 
mal, while the singular values are non-negative and or¬ 
dered in size, Ai > A 2 > ... > Aat > 0. 

Such a singular value decomposition can be applied to 
any matrix M, say of rank TV, where the size of the N 
non-vanishing singular values Ai signifies the significance 
of a singular vector. Conversely, given the decomposition 
of the matrix, it can be approximated by a matrix of rank 
X by only keeping the y largest singular values. Actually 
this matrix is the best approximation of M by a matrix 
of rank y with respect to the least square error. Whether 
this is a good approximation can be directly estimated 
from the size of the discarded singular values with respect 
to the kept ones, in particular the largest. 

In the context of tensor network renormalization, the 
matrices U and V are variable transformations, in fact 
isometries. They map the fine degrees of freedom, en¬ 
coded in the indices a,6, etc., into coarse degrees of free- 


index, such that the entries of the tensor are arranged as the 
indices of a matrix. 


dom that are (for the moment) labelled by the number 
of the singular value. They also give the intermediate 
3-valent tensors: 

. (B4) 

Truncating the number of singular values results in a 
truncation on the index range on the new edges of the in¬ 
termediate tensor network. As a final step, the fine edges, 
labelled by the original indices a, b ,..., are contracted and 
one defines a new tensor T': 

^a'b'c'd' ^ ^ {,^l)dc,c'{^2^ba,a'iS3^ad,c'{.^4)cb,a' 

a,b,c,d 

(B5) 

The new tensor network consisting of effective tensors 
T' is coarser and tilted by 45 degrees. From here on 
the procedure is iterated; the system has reached a fixed 
point if the tensor does not change under consecutive 
iterations, described by a fixed point tensor T*. 

While this truncation solves the practical problem of 
exponentially growing index ranges, the interpretation of 
the new degrees of freedom is unclear. In each itera¬ 
tion of the algorithm, the fine variables get redefined and 
lose their meaning (in terms of the original variables) if 
one forgets these variable redefinitions, yet even keeping 
them is not very feasible after many iterations. More¬ 
over, if a system possesses symmetries, e.g. from an un¬ 
derlying (quantum) group, these symmetries will also be 
present in the tensor, e.g. the Z 2 delta-function present 
on each vertex in the Ising model. This symmetry will 
also survive the coarse graining procedure, but may not 
be obvious to identify after several variable redefinitions. 
Hence it is preferable to use an algorithm that explicitly 
preserves the symmetries. 

On the level of the matrix that is to be split one 
can achieve this by considering representation theory of 
the underlying group(s). Since we are considering in¬ 
tertwiners, i.e. elements of the invariant subspace of a 
tensor product of representation vector spaces, on the 
vertices, the irreducible representations on the adjacent 
edges must couple to the trivial representation, e.g. k = 0 
for Z 2 . This 4-valent intertwiner space on the vertices 
can be split into 3-valent intertwiner spaces, analogous 
to the splitting of the tensor network. The new edge 
connecting the 3-valent vertices carries a new represen¬ 
tation, which for Z 2 is uniquely determined. If the 4- 
valent intertwiner space is not one dimensional, as for 
SU(2)fe, the intermediate label can take multiple values. 
Conversely, a specific intermediate label can allow sev¬ 
eral pairs of representations on the fine edges. Crucially, 
the splitted intertwiner admits the same configurations 
of exterior edges as the original one, but allows them to 
be arranged according to the intermediate label. 

The same holds for the matrices M subject to the SVD. 
Their entries can be rearranged according to the inter¬ 
mediate label, which turns them into a block diagonal 
form. It is thus sufhcient to only consider the blocks for 
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each intermediate label and neglect the forbidden conhg- 
urations. This not only has computational advantages, 
since we can perform more SVDs on smaller matrices, but 
also interpretative ones. We call these blocks intertwiner 
channels, see also [53, 54]: Their label will be endowed 
onto the new edge of the coarser tensor network, which 
is the same label as the original model. Thus we also ex¬ 
plicitly preserve the intertwiner structure on the vertices. 
The algorithm also slightly changes because we perform 
a single SVD for each block. In order to obtain a correct 
truncation of the x largest singular values, it is neces¬ 
sary to compare all values from all blocks and take the x 
largest of them. As a result, some intertwiner labels can 
appear with a multiplicity larger than one and lead to a 
generalized class of models with respect to the original 
one. 

In the system at hand this block diagonal form de¬ 
pends on two parameters, a SU(2)fe representation ji and 
a Z 2 representation ki. For 811(2)^ this form is computed 
by expanding it in a particular recoupling basis, essen¬ 
tially of Clebsch-Gordan coefficients, which capture the 


dependence on the magnetic indices which is unchanged 
during the coarse graining procedure. This idea has been 
developed for (finite) non-Abelian groups in [53] and has 
been derived in great detail for quantum groups in [54] 
and will not be repeated here. We add a small section 
in the appendix A to provide a short tour for the inter¬ 
ested reader. For the Ising model, it is straightforward 
and can be found in detail in [51, 52], which deals with 
finite Abelian groups Zg with q > 2. Here we will briefly 
present the idea: 

As discussed above, the Z 2 representations ke meet¬ 
ing at an edge have to sum to zero (modulo 2 ), which 
is encoded in the Z 2 delta-function on the vertex. If 
we split the 4-valent vertex into two 3-valent ones, we 
replace one delta-function by two and introduce an ad¬ 
ditional variable ki. The delta-functions enforce that 
ki + k 2 = ki = k^ + k 4 , such that solving one constraint 
for ki restores the original delta-function. Conversely, if 
both (5-functions are solved for e.g. k 2 and ^ 4 , the tensor 
is in block diagonal form with parameter ki. 

To sum up this discussion, the tensors T that are meant 
to be split, can be written in the following form: 


J2 3i 


{mj, {fcj) = {h}) (5(2)(fci + k2- fcs) <5(2) (fcg + k4- fcs) 


J5 


J3 Ji 


32 3l 

Tm, {mj, {fcj) = {fcj) + h- fee) 5(2) (fc4 + fci - fee) 


33 3i 


(B 6 ) 


(B7) 


The two different basis correspond to the two different 
splittings of the tensors, where js, k^ and Jq, fcg label 
the intermediate labels respectively, the graphs pictori- 
ally encode the dependence on magnetic indices. In prin¬ 
ciple, the delta-functions encoding the constraints of the 


Ising model can be implicitly included in the T^, yet we 
write them out to underline the fact that this symmetry 
is preserved by the algorithm. In the next step, the SVD 
is applied to T for each choice of intermediate labels. One 
computes the following 3-valent tensors: 
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J2 Jl 
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45 



(B9) 


iS3Yr^f^i{I}{2,3},l) + h- h) 

h 

jl 

J6 ^ 
44 


J3 44 
42 


(^4)^^^({/}ri},*) +h- k,) ■ 


(BIO) 


(Bll) 


The I in the argument of Si, i = 1 ,..., 4, summarizes both SU(2)jt and Z 2 representations. Note that both the Clebsch- 
Gordan coefficients and the Z 2 delta-functions split trivially due to our choice of basis. The single non-trivial step 
is the SVD on the T^. As the final step, these four 3-valent tensor are combined according to the fig. 8 to form the 
new tensor T'. To obtain the recurrence relation for T', we have to contract the magnetic indices and obtain (see also 
appendix A): 


a,b,c,d {m} 


4241 
43 44 


45 


X (^2)^^^(W{M}>*l) (^4)^^^({/}{6 .c}G2) 

X (^l)^4"MU}{rf.e},*3)(^3)^^f ({/}{a.d}A4) • (B12) 

After substituting the equations for all Si and contracting all magnetic indices, we obtain the final recurrence relation: 

(_l)4c+4a+45 
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where we have used the following identity of Clebsch-Gordan coefficients that is proven in [54] and appendix A: 


42 


43 



(- 1 )^' 


4c+4a+45 (_X)4l+42+43+i4 


djb i/djbdj~d 


Jc 3a 45 
4l 42 jb 


3c 3a 45 
44 43 3d 
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This concludes the discussion of this particular tensor 


network algorithm. Interestingly, it can be simplified fur- 
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ther by ‘cutting it in half’ and using instead an algorithm 
based on 3-valent vertices, which is straightforward to 


derive from the 4-valent one and the subject of section 
III. 
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